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Abstract 


Superior  performance  is  a  key  goal  in  weapons  design.  Missiles  must  fly  faster  and  farther,  and  be 
more  agile,  compact  and  lightweight.  The  United  States  have  been  very  active  in  the  development 
of  smart  structures  for  military  applications  such  as  missile  guidance.  For  example,  DARPA  and 
the  Air  Force  Office  of  Scientific  Research  are  sponsoring  the  development  of  miniaturized  active 
flow  technologies  to  achieve  far-term  objectives  for  smart  bombs  and  missiles.  Their  principal 
concepts  for  altering  the  flow  around  a  body  are  centred  on  miniaturized  devices  that  are 
embedded  in  the  missile  skin  and/or  airframe.  The  objective  of  this  project  is  to  demonstrate  with 
hardware  and  software:  1)  engineered  micro-flow  effectors  of  specific  geometry  and  placement 
that  are  able  to  produce  controllable  aerodynamic  forces  on  a  missile  under  supersonic  conditions 
or  a  delta  wing  under  subsonic  conditions,  2)  microactuators  that  are  able  to  meet  the  force, 
kinematic  and  thermal  requirements  set  by  the  aerothermal  environment,  and  3)  control 
algorithms  running  on  non-flightweight  electronics  and  feedback  sensors  that  take  in  command 
signals  and  output  appropriate  actuator  drive  signals  to  produce  the  desired  aerodynamic  force  on 
a  wind  or  water  tunnel  delta  wing  model.  This  document  records  the  progress  made  by  the  project 
members  for  fiscal  year  2005/2006  in  the  areas  of  missile  aerodynamics,  delta  wing 
aerodynamics,  microactuator  modeling,  control  synthesis  and  micro-fabrication. 


Resume 


Une  performance  superieure  est  un  objectif  fondamental  dans  la  conception  d'armes.  Les  missiles 
doivent  voler  plus  rapidement,  sur  de  plus  grandes  distances  et  etre  plus  agiles,  plus  compacts  et 
plus  legers.  Les  Etats-Unis  ont  ete  tres  actifs  dans  le  developpement  de  structures  intelligentes 
pour  des  applications  militaires  telles  que  le  guidage  des  missiles.  Par  exemple,  DARPA  et  YAir 
Force  Office  of  Scientific  Research  commanditent  le  developpement  de  technologies 
d’ecoulement  actif  en  vue  d’objectifs  a  long  terme  pour  les  bombes  et  missiles  intelligents.  Leurs 
principaux  concepts  pour  modifier  fecoulement  autour  d'un  corps  portent  sur  les  dispositifs 
miniaturises  qui  font  partie  du  revetement  et/ou  de  la  cellule  du  missile.  L'objectif  de  ce  projet  est 
de  demontrer  :  1)  que  les  microsurfaces  d’une  forme  particuliere  et  leur  mise  en  place  peuvent 
produire  des  forces  aerodynamiques  controlables  sur  un  missile  dans  des  conditions 
supersoniques  ou  une  aile  delta  dans  des  conditions  subsoniques;  2)  que  les  microactionneurs 
peuvent  repondre  aux  exigences  de  force,  cinematique  et  thermique;  et  3)  que  les  algorithmes  de 
controle  fonctionnant  sur  les  equipements  du  laboratoire  captent  les  signaux  d’ entree  et  generent 
des  signaux  qui  produisent  la  force  aerodynamique  desiree  sur  un  modele  d'aile  delta  dans  un 
tunnel  aerodynamique  ou  hydrodynamique.  Ce  document  fait  etat  des  progres  accomplis  par  les 
participants  au  projet  durant  l’annee  financiere  2005-2006  dans  les  domaines  de  faerodynamique 
du  missile,  de  fade  delta,  de  la  modelisation  du  microactionneur,  de  la  synthese  du  controle  et  de 
la  microfabrication. 


DRDC  Valcartier  TN  2006-178 


This  page  intentionally  left  blank. 


ii 


DRDC  Valcartier  TN  2006-178 


Table  of  Contents 


Abstract/Resume . i 

Table  of  Contents . iii 

1.  Introduction . 1-1 

2.  Lessons  Learned  about  Multi-disciplinary  Interaction . 2-1 

2.1  Participant  Impressions . 2-2 

3.  Cumulative  Action  Items . 3-1 

4.  Missile  and  Delta  Wing  Aerodynamics . 4.1-0 

4.1  CFD  and  Experimental  Investigation  on  Side-Force  Generation . 4.1-0 

4.2  A  CFD  Parametric  Study  on  Missile  Yawing  Control  using  Nose-Mounted  Flow 

Effectors . 4.2-0 

4.3  Impact  of  Nose-mounted  Micro-structures  on  the  Aerodynamics  of  a  Generic 

Finned  Missile . 4.3-0 

4.4  Mini-scale  PIV  Measurements  during  Unsteady  Excitation  of  Mini-Leading  Edge 

Flaps  Used  for  Separation  Control . 4.4-0 

5.  Microactuator  Development . 5.1-0 

5 . 1  SMA  Antagonistic  Actuator  Modeling  and  Validation  Results . 5.1-0 

5.2  Compliant  Flow  Effector  Design  for  Missile  Side  Force  Control . 5.2-0 

6.  Micro-Fabrication  Technology . 6.1-0 

6.1  Precision  Fabrication  of  Active  Flow  Effectors  for  Supersonic  Missile  Control 
. 6.1-0 

7.  Control  Synthesis . 7.1-0 

7. 1  Position  Control  of  SMA  Actuation . 7.1-0 

7.2  Nonlinear  Control  Methods  for  Delta  Wing  Systems  with  SMA  Actuator 

Dynamics . 7.2-0 

8.  Participant  Information . 8-1 

9.  Bibliography . 9-1 

Distribution  List 

DRDC  Valcartier  TN  2006-1 78  iii 


This  page  intentionally  left  blank. 


IV 


DRDC  Valcartier  TN  2006-178 


1.  Introduction 


Superior  performance  is  a  key  goal  in  weapons  design.  Missiles  must  fly  faster  and  farther,  and  be 
more  agile,  compact  and  lightweight.  The  United  States  has  been  very  active  in  the  development 
of  smart  structures  for  military  applications  such  as  missile  guidance.  For  example,  DARPA  and 
the  Air  Force  Office  of  Scientific  Research  are  sponsoring  the  development  of  miniaturized  active 
flow  technologies  to  achieve  far-term  objectives  for  smart  bombs  and  missiles.  Their  principal 
concepts  for  altering  the  flow  around  a  body  are  centred  on  miniaturized  devices  such  as 
magnetic  micro  flap  actuators,  bubble  actuators,  micro-synthetic  jets,  micro-jet  actuators  and  Co- 
Located  Actuator  and  Sensor  (CL AS)  units. 

The  success  of  using  smart  structures  for  active  flow  control  depends  on  a  micro-flow  effector’s 
ability  to  influence  the  macroscopic  flow  around  the  missile.  Missiles  with  slender  forebodies 
face  significant  yawing  moments  under  asymmetric  vortices.  Controlling  forebody  vortex 
asymmetry  is  dependent  on  the  sensitivity  of  the  asymmetric  flow  to  the  distance  between  the 
micro-flow  effector  and  the  nose  tip.  The  closer  the  micro-flow  effector  is  to  the  nose  tip,  the 
lower  the  power  required  to  trigger  flow  changes.  Several  techniques  for  manipulating  forebody 
vortices  have  been  investigated.  Most  methods  are  essentially  steady  schemes  that  produce  quasi¬ 
steady  loads  by  forcing  the  forebody  vortices  into  desired  positions. 

The  Defence  R&D  Canada  project  “Supersonic  Missile  Flight  Control  by  Manipulation  of  the 
Flow  Structure  using  Micro-Actuated  Surfaces”  was  officially  started  on  1  April  2003.  The 
project  termination  date  is  31  March  2006.  The  objective  of  the  project  is  to  demonstrate  with 
hardware  and  software  the  following  items: 

a)  Engineered  micro-flow  effectors  of  specific  geometry  and  placement  that  are  able  to 
produce  controllable  aerodynamic  forces  on  a  missile  under  supersonic  conditions  or  a  delta 
wing  under  subsonic  conditions. 

b)  Microactuators  that  are  able  to  meet  the  force,  kinematic  and  thermal  requirements  set  by 
the  aerothermal  environment.  The  availability  of  microactuators  necessarily  means  that 
material  and  fabrication  issues  have  been  resolved  to  the  extent  that  prototypes  can  be  made. 

c)  Control  algorithms  running  on  non-flightweight  electronics  and  feedback  sensors  that 
take  in  command  signals  and  output  appropriate  actuator  drive  signals  to  produce  the  desired 
aerodynamic  force  on  a  wind  or  water  tunnel  delta  wing  model. 

A  meeting  was  held  on  19  April  2006  at  DRDC  Valcartier  to  make  formal  presentations  that 
informed  the  participants  on  the  progress  that  was  made  from  1  April  2005  to  31  March  2006.  As 
a  secondary  objective,  the  meeting  was  held  to  promote  a  better  understanding  on  the  technical 
issues  that  must  be  addressed  in  future  work  as  the  team  works  towards  achieving  the  project 
objectives. 

All  participants  were  asked  to  provide  a  250  word  abstract  and  a  PowerPoint  presentation 
containing  the  following  elements  as  they  were  applicable: 

1)  title  page  stating  title  of  presentation,  contributors,  name  of  organization  and  date 
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2)  previous  studies  in  the  open  literature, 

3)  geometry  studied, 

4)  test  set  up  and  instrumentation, 

5)  modeling  assumptions, 

6)  governing  equations  of  phenomenon, 

7)  parametric  study, 

8)  major  observations  and  problems, 

9)  conclusions, 

10)  future  work. 

This  documents  collects  the  presentations  made  by  the  project  members  in  the  areas  of  missile 
aerodynamics,  delta  wing  aerodynamics,  microactuator  modeling,  control  synthesis  and  micro¬ 
fabrication.  The  abstracts  provided  by  the  participants  precede  their  presentations. 
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2.  Lessons  Learned  about  Multi-disciplinary 
Interaction 


Development  of  active  flow  control  technology  requires  advances  in  the  areas  of  aerodynamics, 
microactuator  design  and  fabrication  and  control  synthesis  to  construct  a  workable  active  flow 
control  system.  The  hierarchy  shown  in  Figure  1  described  best  the  interrelationships  and  flow  of 
the  development  process.  Because  active  flow  attempts  to  alter  the  flow  structure,  the 
fundamental  problem  in  the  case  of  aerovehicles  is  one  of  aerodynamics.  The  success  of  active 
flow  control  technology  rests  on  the  proper  understanding  of  the  interaction  between  parameters 


Aerodynamics 


Figure  1.  Information  flow  between  technologies 

that  affect  the  nonlinear  flow  phenomenon  and  the  resultant  aerodynamic  forces.  The 
microactuator  design  can  proceed  once  the  aerodynamic  loads  and  the  required  deflections  and 
actuation  frequencies  are  known.  Actuator  modeling  activities  are  centred  on  materials  selection, 
characterization  and  mechanism  synthesis.  Control  synthesis  necessarily  requires  a  mathematical 
description  of  the  processes  that  it  aims  to  control.  In  this  case,  control  algorithm  development 
aims  to  optimize  a  highly  nonlinear,  time-dependent  airflow-microactuator  response  to  a  given 
command  input.  The  system  model  is  used  to  evaluate  the  performance  of  a  proposed  aero- 
mechanical-control  system  and  can  provide  insight  as  to  the  areas  where  performance  gains  can 
be  obtained.  The  micro-fabrication  activity  is  dependent  on  the  progress  made  in  the 
microactuator  design.  With  the  great  strides  made  in  the  fabrication  of  micro-electromechanical 
systems,  the  problem  for  active  flow  control  devices  is  centred  on  the  development  of  machining 
and  joining  techniques  to  make  a  miniature  multi-material,  multi-component  device. 

The  progress  in  the  second  year  of  the  project  has  been  steady  though  arduous.  Each  technology 
group  worked  in  semi-isolation  in  order  to  concentrate  on  the  detailed  model  development  or 
experimental  measurement  particular  to  the  technology.  This  step  was  necessary  so  that  there  is  a 
high  confidence  in  the  data  and  tools  that  will  be  used  in  the  system  level  trade-off  studies. 


DRDC  Valcartier  TN  2006-178 


2-1 


In  the  final  year,  the  experimental  missile  aerodynamic  work  was  reduced  to  allow  the  CFD 
studies  time  to  examine  the  flow  physics  in  detail.  With  the  position  and  geometry  of  the  flow 
effectors  determined,  it  was  possible  to  design  a  dynamic  flow  effector  because  the  capability  of 
the  SMA  actuator  was  also  known.  The  process  followed  a  linear  progression  moving  from  flow 
effector  design  to  micro-fabrication  to  sub-scale  assembly  and  then  to  control  synthesis.  The 
preliminary  studies  carried  out  for  micro-fabrication  and  control  algorithm  development 
permitted  insights  to  be  made  into  the  factors  that  affected  the  performance  in  each  of  the 
technology  areas. 

2.1  Participant  Impressions 


At  the  end  of  the  2005/2006  review  meeting,  a  questionnaire  was  given  to  the  participants  to 
record  their  impressions  on  lessons  learned,  group  interaction  and  ideas  for  potential  follow-on 
projects.  The  following  paragraphs  summarize  their  assessments. 

1.  What  were  the  ‘lessons  learned’  for  your  particular  technology? 

The  flow  effector  position  can  influence  the  vortex  system  so  differently.  Excellent  project  to 
have  multiple  disciplines  together  to  understand  the  problem/project  in  a  comprehensive 
way. 

CFD  results  are  dependent  on  the  turbulence  model  used.  A  more  exhaustive  study  using 
different  turbulence  models  and  CFD  codes  would  provide  a  greater  level  of  confidence  in 
the  results. 

You  never  know  until  you  try  but  all  reasonable  results  come  from  careful  thinking  and  deep 
understanding  of  the  principles  of  flow  behaviour. 

From  the  aerodynamic  point  of  view,  the  application  of  the  micro-flow  effector  to  a 
supersonic  missile  appears  limited.  Vortices  can  only  be  affected  at  low  speeds  or  low 
supersonic  Mach  number.  Vortices  form  at  angles  of  attack  between  10  and  20  deg.  It  is  not 
obvious  to  have  a  missile  flying  at  such  an  attitude.  Micro-flow  effectors  might  have  more 
application  for  cancelling  side  forces  than  to  generate  side  forces. 

For  SMA’s,  cooling  time  constants  of  the  material  limit  the  actuation  frequency.  Forced 
cooling  can  be  used  however  the  actuator  then  loses  its  simplicity  of  design. 

To  achieve  a  high  actuation  rate,  the  configuration  of  the  actuator  must  be  considered  in 
addition  to  actuator  material.  The  antagonistic  arrangement  is  a  good  example.  Maybe 
another  arrangement  could  be  better. 

Experimentation  for  SMA  control  was  the  main  issue.  Several  iterations  between  theory  and 
practice  were  and  are  still  necessary. 

Require  experimental  results  for  control  proofs.  Sophisticated  control  methods  are  needed  to 
achieve  satisfactory  performance.  Delta  wing  problem  is  a  good  demonstration  of  control 
synthesis  development. 
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Experimentation  is  very  important  for  active  flow  control  development. 

Timely  development  of  micro-milling  technology  allowed  fabrication  of  flow  effectors  with 
desired  accuracy.  For  higher  production  rates,  laser  and  milling  technique  would  work 
effectively.  For  features  sizes  less  than  100  micron  and  acute  angles,  lasers  may  still  be  an 
effective  tool. 

2.  What  were  the  ‘lessons  learned’  while  working  in  a  multi-disciplinary  group? 

Could  start  with  CFD  conceptual  computations  then  design  experiments  for  several  typical 
cases.  Detailed  CFD  could  be  re-run  afterwards  to  corroborate  results. 

The  flow  of  information  begins  with  the  aerodynamic  results.  A  viable  aerodynamic  control 
method  is  required  to  set  the  context  for  the  structures,  controls  and  fabrication  efforts. 

Two  heads  are  better  than  one.  It  is  very  important  to  have  a  multi-disciplinary  team  in  one 
group  to  make  a  project  work  and  to  ensure  that  people  are  benefltting  from  each  other’s 
expertise. 

It  is  important  to  start  the  integration  of  the  various  technologies  early  on  in  the  project  in 
order  to  have  enough  time  to  face  the  inherent  problems  associated  with  the  various 
technologies. 

Controlling  the  actuator  depends  on  the  mechanical/electrical  design,  the  aerodynamics  and 
the  control  algorithms.  Integrated  design  is  necessary  to  achieve  good  performance  of  multi- 
component  systems. 

It  is  difficult  to  work  on  a  problem  with  variables  that  depend  on  other  technology  domains 
especially  when  those  variables  are  essential  to  advance  our  own  work. 

It  was  interesting  to  see  how  other  disciplines  approached  technical  problems. 

Early  development  of  the  delta  wing  vortex  breakdown  model  was  very  useful  for  the 
development  of  roll  control  algorithms. 

Difficult  to  understand  other  technologies.  Not  everyone  understood  the  need  for  dynamic 
feedback  control.  Sometimes,  each  technology  assigned  different  meanings  to  similar  words. 

Understanding  of  the  aerodynamics,  modeling  the  flow  patterns  and  the  design  and  control  of 
the  actuators  laid  out  a  clear  path  towards  product  development  although  there  is  still 
significant  practical  work  that  remains  to  be  done. 

3.  What  could  the  group  have  done  better  during  the  past  three  years  to  advance  the 
technology  beyond  the  point  that  we  are  at  today? 

Some  flow  visualization  studies  in  addition  to  CFD  analyses  would  have  helped  understand 
the  flow  mechanism  better. 
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The  decision  to  carry  out  CFD  studies  on  the  missile  and  delta  wing  configuration  allowed  us 
to  minimize  the  risks  associated  with  the  project.  However,  a  review  after  the  first  year’s 
effort  to  decide  which  platform  to  focus  on  could  have  permitted  greater  comprehension  of 
the  flow  physics  on  the  chosen  platform. 

Have  group  review  meetings  twice  a  year  instead  of  just  one  meeting. 

A  better  literature  review  on  missile  aerodynamics  could  have  allowed  a  better  application  of 
the  micro-flow  effector  to  be  developed. 

A  common  analysis  model  covering  the  aerodynamic-flow  effector  interaction,  the  compliant 
mechanism-SMA  dynamics  and  the  control  algorithm  synthesis  would  have  permitted  more 
concepts  and  systems  issues  to  be  examined. 

The  group  worked  well  together  during  the  last  three  years  to  achieve  the  objectives. 
However,  the  time  required  to  complete  all  objectives  may  have  been  underestimated  which 
would  explain  why  the  project  could  not  be  completely  finished  by  3 1  March  2006. 

Having  monthly  review  meetings  is  required  to  stay  on  schedule. 

More  interaction  between  the  design  of  the  compliant  mechanism  and  the  control  group 
would  have  been  helpful. 

More  understanding  of  other  people’s  technologies.  A  more  focused  approach  with  fewer 
topics  would  have  helped. 

A  student  internship  program  could  have  been  helpful. 

Early  input  on  aerodynamics  and  modeling  studies  would  have  taken  the  project  further  in  the 
completion  of  the  first  prototype  integrated  device.  Most  the  results  achieved  remain  at  the 
fundamental  understanding/scientific  level. 

4.  What  follow-on  studies  can  you  forsee  as  a  result  of  the  technology  developed  in  this 
project? 

Continue  to  examine  the  flow  effector  and  missile  fin  configuration  to  examine  the 
interaction  between  the  fin  and  the  vortex  caused  by  the  flow  effector.  Look  at  Reynolds  and 
Mach  number  effects.  Continue  CFD  study  of  delta  wing. 

Control  of  the  boundary  layer  on  flapping  wings  or  delta  wings. 

Drag  reduction  of  projectiles  at  high  angles  of  attack. 

Application  of  micro-flow  effectors  to  low  velocity  ammunition. 

Long  range  unspun  artillery  projectile  control. 
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Actuators  for  munitions  must  remain  compact  yet  be  able  to  produce  displacements  in  the 
millimeter  range  at  frequencies  up  to  100  Hz.  Piezoceramic  materials  coupled  to  a  compliant 
mechanism  is  an  area  that  could  be  investigated. 

Investigate  MEMS  to  speed  up  dynamic  response. 

Solenoid  actuators  for  flow  control  and  further  experimentation  on  the  delta  wing. 

Outer  loop  flow  control  integrated  with  NRC  delta  wing  testbed. 

Examine  implementation  issues  such  as  real-time  performance  of  receding  horizon  control  for 
outer  loop  control  of  delta  wing. 

Would  like  to  see  the  validation  of  the  prototype  missile  flow  effector  device. 


Figure  2.  Project  participants  (r  to  I)  Olivier  Boissonneault,  Brandon  Gordon,  Suwas  Nikumb, 
Daniel  Corriveau,  Nicolas  Lechevin,  Camille  Alain  Rabbath,  Xingzhong  Huang,  Suzhen  Chen, 

Nicolas  Hamel,  Frank  Wong 
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3.  Cumulative  Action  Items 


Legend  for  groups: 

DR  DC  -aero:  Daniel  Corriveau,  Nicolas  Hamel,  Francois  Lesage 

DRDC-struct:  Frank  Wong,  Olivier  Boissonneault  (ETS) 

Numerica:  Pierre  Antoine  Ranville,  Pierre  Gosselin 

IAR-cfd:  Suzhen  Chen,  Stuart  Mcllwain,  Mahmood  Khalid 

IMTI:  Suwas  Nikumb,  Dan  Zhang 

ETS:  Patrick  Terriault 

Concordia:  Wei  Wei  Liu,  Brandon  Gordon 

Legend  for  individuals: 

IAR-xzh:  Xing  Zhong  Huang 
IAR  -sc:  Suzhen  Chen 
IMTI-sn:  Suwas  Nikumb 
DRDC-dc:  Daniel  Corriveau 
DRDC-car:  Camille  Alain  Rabbath 
DRDC-fw:  Frank  Wong 


April  2005  Action  Items 

Y 0405-1 .  Continue  with  pressure  measurements  and  oil  or  vapor  screen 

visualization  studies  on  the  DRDC-B1AC2R  model  with  and  without  fins.  Determine 
how  separation  line  changes  with  flow  effector  position. 

Action:  DRDC-dc  (May-Sep  2005) 

Y0405-2.  Compare  pressure  measurements  and  flow  visualization  results  with 
current  set  of  CFD  results. 

Action:  DRDC-dc,  DRDC-nh,  IAR-sc  (May- July  2005) 

Y0405-3.  Carry  out  additional  CFD  analysis  for  the  finless  DRDC-B1 AC2R 
between  15  deg.  to  20  deg.  to  determine  where  the  peak  side  force  is  situated.  Compare 
pressure  distribution  with  IAR  results. 

Action:  DRDC-nh  (April-May  2005) 

Y0405-4.  Produce  joint  DRDC-IAR  paper  comparing  CFD  to  wind  tunnel  results. 
Propose  explanation  on  how  the  flow  effector  generates  the  observed  side  forces. 
Action:  DRDC-dc,  IAR-sc 

Y0405-5.  Produce  joint  DRDC-IAR  paper  comparing  results  from  Fluent  and 
SPARC  CFD  codes  and  various  turbulence  models. 


DRDC  Valcartier  TN  2006-178 


3-1 


Action:  DRDC-nh,  IAR-sc 


Y0405-6.  Carry  out  CFD  study  on  DRDC-B1AC2R  with  fins.  Decide  on  analysis 
matrix.  Coordinate  with  wind  tunnel  test  matrix. 

Action:  DRDC-nh 

Y0405-7.  Decide  on  the  best  way  to  visualize  the  separation  line  from  the  CFD 
results. 

Action:  DRDC-nh,  IAR-sc  (May  2005) 

Y0405-8.  Characterize  Dynalloy  SMA.  Construct  small-scale  test  bench. 

Action:  DRDC-ob  (May  2005) 

Y0405-9.  Examine  ability  of  SMA  model  to  predict  response  under  complex  force 
profile. 

Action:  DRDC-ob  (June-July  2005) 

Y0405-10.  Investigate  methods  to  increase  frequency  response  of  the  SMA  wire 
given  force  requirements. 

Action:  DRDC-ob,  DRDC-fw  (Aug-Sep  2005) 

Y0405-1 1 .  Examine  impact  on  the  control  strategy  for  a  semi-active  actuator. 
Aeroloading  is  assumed  negligible.  F.  Wong  will  provide  spring  constant. 

Action:  DRDC-nl,  DRDC-car  (May  2005) 

Y0405-12.  Clarify  what  aerodynamic  data  was  provided  to  Concordia  and  agree  on 
what  additional  data  will  be  provided. 

Action:  Concordia,  IAR-xzh 

Y0405-13.  Develop  method  to  analyze  compliant  beam.  Integrate  method  with  SMA 
model  to  determine  appropriate  size  for  beam  elements. 

Action:  DRDC-fw,  DRDC-ob  (May- June  2005) 

Y 0405-14.  Develop  manufacturing  plan  for  micro-flow  effector. 

Action:  DRDC-fw,  IMTI-sn  (June  2005) 

Y0405-15.  Preliminary  test  to  integrate  MatLab  control  law  in  LabView. 

Action:  DRDC-nl,  DRDC-car  (May- June  2005) 

Y0405-16.  Integrate  control  law  in  macro-  and  micro-test  bench  environment. 
Explore  ability  to  track  command  signal. 

Action:  DRDC-car,  DRDC-nl,  DRDC-fw,  DRDC-ob  (July  2005) 

Y0405-17.  Micro-PIV  study  on  2D  wing  with  pressure  taps. 
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Action:  IAR-xzh  (May  2005) 


Y0405-18.  Provide  delta  wing  drawings  to  DRDC.  DRDC  will  provide  minimum 
wing  thickness  based  on  actuator  requirements. 

Action:  IAR-xzh,  DRDC-fw  (May  2005) 

Y0405-19.  Coordinate  fabrication  of  delta  wing  between  IAR  and  IMTI. 

Action:  DRDC-fw 


Y0405-20.  Investigate  possibly  of  loaning  U.  Laval  micro-PIV  system  to  IAR  to 
continue  delta  wing  study. 

Action:  DRDC-nh  (April  2005) 


Y 0405-2 1 .  Prepare  validation  test  plan  for  missile  model  with  micro-actuated  flow 

effectors. 

Action:  DRDC 

Y0405-22.  Prepare  validation  test  plan  for  delta  wing  model  with  micro-actuated 
leading  edge  fence. 

Action:  DRDC-fw,  DRDC-car,  IAR-xzh 


Y0405-23.  Provide  list  of  papers  that  will  be  produced  for  2005/2006  including 
conference  details. 

Action:  All 


April  2004  Action  Items 

Y 0304- 1 .  Advance  individual  technology  development  from  April  2004  to  February 

2005  to  the  point  that  integration  issues  can  be  addressed  in  March  2005. 

Action:  All 

Y0304-2.  Missile  geometry  will  be  referred  to  as  the  DRDC-B1AC2R  geometry. 
Flow  alteration  devices  will  generally  be  referred  to  as  flow  effectors. 

Action:  All 

Y0304-3.  A  common  definition  for  the  missile  coordinate  system  will  be  agreed 
upon. 

Action:  DRDC-aero,  Numerica,  IAR-cfd 

Y0304-4.  A  standard,  compact  nomenclature  to  describe  missile  flow  effector 
number,  placement  and  orientation  will  be  developed  for  use  by  all  participants. 

Action:  DRDC-dc 
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Y0304-5.  DRDC  will  continue  using  blunt-edged  rectangular  flow  effectors  in  its 
wind  tunnel  tests  unless  CFD  studies  show  there  is  an  advantage  to  using  a  knife-edged 
flow  effector. 

Action:  DRDC-dc 

Y0304-6.  Verify  Numerica  results  and  the  impact  of  knife-edged  flow  effectors. 
Action:  DRDC-nh,  Numerica,  IAR-cfd 

Y0304-6. 1  Common  conditions:  15  deg  AOA,  flow  effector  placed  at  front 
station,  no  base  drag,  ffeestream  conditions  M  =  1.5,  26kPa,  206K 

Y0304-6.2  DRDC  will  perform  a  grid  refinement  study  of  the  short  conical 
nose  body  with  one  flow  effector  at  45  deg  orientation  to  verify  if  the  Numerica 
CFD  result  is  grid  independent. 

Y0304-6.3  Numerica  will  modify  the  CFD  mesh  for  the  short  conical  nose 
body  with  one  flow  effector  at  45  deg  orientation  to  include  a  knife-edge  on  the 
leading  and  trailing  edges  of  the  flow  effector.  Each  knife-edge  will  occupy  one- 
third  of  the  flow  effector  chord. 

Y0304-6.4  IAR  will  verify  the  Numerica  blunt-edged  flow  effector  result  by 
analyzing  a  short  conical  nose  body  with  one  blunt-edged  rectangular  flow  effector 
oriented  at  45  deg. 

Y0304-7.  If  CFD  studies  show  there  is  no  difference  between  a  blunt-edged  and 
knife-edge  flow  effector,  decide  whether  it  is  acceptable  to  compare  CFD  results  derived 
from  a  knife-edged  flow  effector  with  wind  tunnel  data  derived  from  a  blunt-edged  flow 
effector. 

Action:  DRDC-aero,  IAR-cfd 

Y0304-8.  Determine  through  wind  tunnel  and  CFD  analyses  on  DRDC-B1AC2R 
with  flow  effectors  whether  it  is  separation  or  vortex  phenomenon  or  a  combination  of 
the  two  that  produces  the  observed  side  forces. 

Action:  DRDC-aero,  Numerica  and  IAR-cfd 
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Y0304-9.  CFD  results  will  be  produced  for  the  test  cases  shown  in  Table  1.  The 
results  will  be  compared  and  conclusions  on  their  consistency  and  trends  will  be  made. 
We  wish  to  gain  confidence  that  the  CFD  analyses  are  able  to  predict  the  observed  wind 
tunnel  trends  so  that  missile  flow  effector  specifications  can  be  generated  for 
microactuator  development. 

Action:  DRDC-aero,  Numerica  and  IAR-cfd 


Table  1  -  CFD  analysis  matrix 

Nose 

short  conical 

Body 

finless  cylindrical 

Flow  Effectors: 

number,  placement,  orientation 

0  (baseline) 

1  front  at  0  deg 

1  front  at  45  deg 

3  front  at  0  deg  with  30  deg  separation 

Angle  of  attack 

0,  10,  15,  20  deg 

Freestream  conditions 

M  =  1.5,  26kPa,  206K 

Required  results 

integrated  forces 
streamlines 
pressure  contours 

Cp  profiles  at  DRDC  specified  x/L  intervals 

velocity  profiles 

skin  friction  for  whole  body 

Y0304-10.  Complete  SMA  linear  actuator  model  and  assemble  test  bench  for 
validation  and  closed-loop  control  tests. 

Action:  DRDC-struct,  DRDC-car 

Y0304-1 1.  Reformulate  SMA  linear  actuator  model  into  a  form  suitable  for  control 
synthesis  development. 

Action:  DRDC-struct,  DRDC-car 

Y03 04- 1 2 .  Discuss  collaboration  on  development  of  compliant  mechanism 

algorithms,  compliant  mechanism  design  and  prototype  fabrication. 

Action:  DRDC-struct,  IMTI,  IAR-sc 

Y0304-13.  Discuss  plan  of  work  for  SMA  machining  and  joining  studies. 

Action:  DRDC-fw,  IMTI-sn,  ETS 

Y0304-14.  IAR  to  provide  load,  deflection,  frequency  requirements  for  delta  wing  to 
DRDC  so  that  DRDC  can  examine  implications  for  SMA  actuator  design. 

Action:  IAR-xzh,  DRDC-struct 
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Y0304-15.  DRDC-aero  to  provide  DRDC-struct  with  preliminary  load  estimate  and 
deflection  requirement  for  missile  flow  effector. 

Action:  DRDC-aero 

Y0304-16.  IAR  will  provide  Cl  database  to  Concordia  for  further  development  of 
delta  wing  plant  model. 

Action:  IAR-xzh,  Concordia 

Y0304-17.  Determine  what  information  aerodynamicists  need  to  provide  for  control 
synthesis  in  a  wind  tunnel  setting. 

Action:  DRDC-aero,  DRDC-car,  DRDC-fw 

Y0304-18.  Determine  what  a  closed  loop  control  of  a  wind  tunnel  missile  equipped 
with  microactuated  flow  effector  experiment  would  look  like. 

Action:  DRDC-aero,  DRDC-fw,  DRDC-car 

Y0304-19.  Prepare  global  timeline  to  coordinate  the  rate  of  progress  in  the  individual 
technologies. 

Action:  DRDC-fw 

Y0304-20.  Examine  possibility  of  performing  trajectory  simulations  using  force 
information  from  microactuated  flow  effectors  to  determine  their  system  effectiveness. 
Action:  DRDC-fw 
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4.  Missile  and  Delta  Wing  Aerodynamics 


4.1  CFD  and  Experimental  Investigation  on  Side-Force  Generation 

presented  by  N.  Hamel,  DRDC  -  Valcartier 

Using  nose-mounted  micro-structures,  it  was  shown  experimentally  and  computationally 
in  earlier  studies  that  a  significant  side  force  could  be  generated  at  moderate  angles-of-attack  and 
in  supersonic  flow  on  a  generic  missile  configuration.  Under  these  flow  conditions,  a  “clean” 
slender  body  would  not  normally  experience  any  lateral  or  side  forces. 

In  these  earlier  investigations,  the  side  force  induced  by  the  presence  of  the  micro¬ 
structures  was  shown  to  peak  for  angles-of-attack  between  15  and  20  degrees.  Furthermore,  it 
was  found  that  the  magnitude  of  the  side  force  generated  was  maximized  by  positioning  the  flow 
effectors  at  an  angular  position  close  to  135  or  225  degrees. 

For  the  current  study,  surface  pressure  measurements  in  the  DRDC  trisonic  wind  tunnel 
as  well  as  additional  numerical  simulations  using  FLUENT®  were  made  in  order  to  determine  the 
mechanism  by  which  the  nose-mounted  micro-structures  contributed  to  the  generation  of  a  side 
force  at  moderate  angles-of-attack  where  otherwise  no  side  force  would  be  present. 

A  vortex  shed  off  from  the  flow  effector  entrains  the  matching  leeside  vortex  to 
eventually  merge  into  a  single  stronger  vortex.  The  newly  formed  vortex,  being  stronger  than  the 
opposite  side  vortex,  is  forced  to  move  upward  above  the  missile’s  leeside.  The  formation  of  the 
flow  effector  vortex  arises  because  of  the  pressure  differential  on  either  side  of  the  flow  effector. 
This  pressure  difference  causes  the  flow  to  wrap  around  the  top  edge  of  the  flow  effector,  thus 
resulting  in  a  vortical  flow. 
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f  Objectives 


Determine  the  mechanism  by  which  the  nose- 
mounted  micro-structures  contributed  to  the  side 
force  generation 

Support  the  wind  tunnel  test  results  and  help  in  the 
interpretation  of  the  results 
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Missile  Configuration 

Baseline  Geometry  With  Conical  Nose 


Aspect  Ratio  L/D  =  13.0 
Nose  Aspect  Ratio  Ln/D  =  3.0 
Conical  Nose 

4  Fins  in  +  Configuration  (Removable) 


All  dimensions  in  mm 
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r\)7  Missile  Configuration 

Flow  Effector  Geometry 


Test  Setup  and  Instrumentation 

DRDC  Trisonic  Wind  Tunnel 


950  m3 
Air  Storage  E 


Choke  Valve 


Aux,  Vacuum  ^  ^ 

Tank  0,6x0 .6  m 

Test  Section 


•  Indraft  wind  tunnel 

•  Test  section:  2’x  2’ 

•  0.2  <  Ma  <  4.0 

•  -20°  <  AOA  <  +20  ° 


»  Test  duration:  5-11  sec 
*  Turnaround  time:  30  min. 
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R07  CFD  Set-up 


Fluent  6. 1 

-  RANS  with  the  coupled  explicit  solver 

-  Realizable  k-epsilon  viscous  model 

-  Enhanced  wall  treatment  with  pressure  option 
Wind  tunnel  freestream  conditions 

-  Mach  Number:  1,5 

-  Static  Pressure:  26  000  Pa  | 


A  &} 

A  w 


'1 
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Surface  Static  Pressure  Coefficient 


No  fms 

Conical  nose  L/D=  3 
1  FE  on  the  front  row 
NED  coordinate  system 

-  Blue  0°  to  90° 

-  Green  90°  to  180° 

-  Red  180°  to  270° 
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r uy  Side  Force  Generation  Mechanism 


•  No  fins 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  225° 

•  Front  row 

•  Ma=  1.5  AoA=  20° 

•  Re/m=  1 5.2x1 06 
•Row  3 
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Side  Force  Generation  Mechanism 


u 

m 


Angular  position  (deg) 


•  No  fins 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  225° 

•  Front  row 

•  Ma=  1.5  AoA=  20° 

•  Re/m=  1 5.2x1 06 
•Middle  of  the  FE 
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ruy  Comparison  of  Results 


Surface  Static  Pressure  Coefficient 


No  fins 

Conical  nose  LN/D=  3.0 
Flow  effectors  centered  at  90° 
Front  row 
Ma=  1.5  AoA=  15° 

Re/m=1 5.2x1 06 
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r^y  Comparison  of  Results 


Surface  Static  Pressure  Coefficient 


flAuulaf  jhm4Ioo 


An®j|3  r  ^  SO 


•  No  fins 

•  Conical 

•  Flow  effectors  centered  at  90° 

•  Front  row 

•  Ma=  1.5  AoA=  15° 

•  Re/m=  1 5.2x1 06 


R  &  D  pour  la  defense  Canada 


Comparison  of  Results 


Surface  Static  Pressure  Coefficient 


•  No  fins 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  90° 

•  Front  row 

•  Ma=  1.5  AoA=  15° 

•  Re/m=  15.2x1 06 
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Comparison  of  Results 

Surface  Static  Pressure  Coefficient 


•  No  fins 

•  Conical 

•  Flow  effectors  centered  at  135° 

•  Front  row 

•  Ma=  1.5  AoA=  15° 

•  Re/m=  1 5.2x1 06 


•V 


Comparison  of  Results 

Surface  Static  Pressure  Coefficient 


06 

f-  _  ** 

OS 

04 

S  02 

E 

■ 

g  g1 

s  o, 

1^1 

a  ° 

■0! 

■Ci 

Row  3 

20  4fl  £fl 

Angular  position  4<f»g] 

•  No  fins 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  135( 

•  Front  row 

•  Ma=  1.5  AoA=  15° 

•  Re/m=  15.2x1 06 


4.1-9 


V 


Comparison  of  Results 

Surface  Static  Pressure  Coefficient 


•  No  fins 

•  Conical 

•  Flow  effectors  centered  at  135° 

•  Front  row 

•  Ma=  1.5  AoA=  15° 

•  Re/m=  1 5.2x1 06 


f\)7  Conclusion 


CFD  is  able  to  match  the  wind  tunnel  pressure 
coefficient  (clean  body) 

CFD  is  able  support  the  wind  tunnel  test  and  help 
in  the  interpretation  of  the  results 

The  nose-mounted  micro-structures  contribution  to 
the  side  force  generation  is  now  understand 
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4.2  A  CFD  Parametric  Study  on  Missile  Yawing  Control  using  Nose- 
Mounted  Flow  Effectors 

presented  by  S.  Chen,  IAR/NRC 

IAR  has  performed  a  CFD  investigation  of  the  aerodynamic  performance  of  blunt-edged 
rectangular  key-shaped  flow  effectors  mounted  at  the  nose  of  a  supersonic  flight  missile.  The 
flow  effectors  (FEs)  were  situated  at  the  front  row  position,  located  at  x/D  =  0.844  from  the  nose 
tip.  The  study  focused  on  four  FE  configurations.  They  were:  a  single  FE  situated  at  =  225° 
(RFE1-225),  a  single  FE  situated  at  =  270°  (RFE1-270),  triple  FEs  spaced  30o  apart  from  the 
middle  one  situated  at  §  =  225°  (RFE3-225),  triple  FEs  situated  at  <|>  =  270°  (RFE3-270). 
Parametric  studies  were  also  performed  considering  the  effects  of  the  locations  (azimuthal  angle 
and  axial  locations),  numbers  and  dimensions  (heights)  of  the  flow  effectors  on  the  side  force 
performance.  The  freestream  flow  conditions  were  Mach  number  Moo  =  1.5,  pressure  Poo  = 

26000  Pa  and  temperature  Too  =  206  K  for  all  simulations.  The  computations  were  run  at  angles 
of  attack  (AOA)  of  0°,  10°,  15°,  17°,  19°  and  20°. 

Multi-block  structured  meshes  were  produced  using  ICEM-CFD  HEXA.  Fine  grids  were 
used  to  resolve  the  viscous  boundary  layer  above  the  missile  surface.  The  average  dimensionless 
normal  distance  of  the  first  grid  point  from  the  wall,  y+,  was  less  than  2.  The  total  number  of  grid 
points  was  about  3  millions.  Steady-flow  numerical  solutions  were  obtained  using  the  CFD  code 
SPARC  with  the  Spalart-Allmaras  one-equation  turbulence  model. 

The  numerical  results  indicated  that  the  maximum  peak  side  forces  occurred  at  AOA  = 

17°  for  all  the  configurations  with  flow  effectors  at  the  front  axial  location  x/D  =  0.844.  Among 
the  configurations  studied,  the  single  FE  situated  at  an  azimuthal  angle  of  225°  was  the  most 
efficient  in  generating  the  side  forces.  The  side  force  gradually  decreased  as  the  FE  position 
moved  to  the  horizontal  symmetry  plane  of  the  body  ((|)  =  270°).  The  front  configuration  was 
more  efficient  than  the  back  configuration  in  influencing  the  vortex  system  behaviour.  The  side 
force  did  increase  with  the  number  of  FEs,  however,  the  advantage  was  not  dramatic.  These 
quantitative  observations  agreed  with  the  wind  tunnel  results  obtained  from  DRDC-Valcartier. 
Regarding  the  quality  of  the  comparison,  the  side  force  coefficients  on  the  starboard  side  of  the 
body  agreed  well  with  the  wind  tunnel  data,  but  those  on  the  port  side  of  the  body,  which  had 
identical  numerical  counterparts,  were  consistently  lower  than  the  wind  tunnel  data.  The 
asymmetrical  features  of  the  wind  tunnel  data  might  have  been  caused  by  the  understandable 
imperfections  in  the  experimental  models  or  the  facilities,  which  could  be  significant  for  the  small 
magnitude  of  the  side  force  coefficients.  The  numerical  normal  and  axial  forces  agreed  with  the 
wind  tunnel  data  well. 

Both  numerical  and  experimental  studies  concluded  that  for  the  studied  cases,  the 
location  of  FEs  was  more  important  than  the  number  of  FEs  in  terms  of  their  contribution  to  the 
side  force.  The  numerical  study  also  showed  that  the  side  forces  generated  by  the  RFE 1-225 
configuration  were  on  average  about  25%  of  the  normal  forces.  The  side  force  gain  was  quite 
remarkable  with  such  a  small  device  installed  at  the  nose  of  the  body.  This  re-confirms  the 
feasibility  of  using  small  devices  at  the  nose  of  the  slender  body  to  control  the  side  forces. 
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Numerical  Descriptions 


Geometry  and  flight  conditions 

•  FEs  located  at  L/D=0.844  (Front  row)  or 
x/D=1.269  (back  row). 

•  Shape:  Rectangular  key 

•  M-1.5,  P=26000Pa,  T=206K, 

RJm  =  1 5.2x1 06. 

AOA  =  0°  to  20 0 


h/D  =0.033 
w/D  =  0.089 
t/D  =0.026 


Numerical  Descriptions 

•  CFD  Solver:  SPARC 

•  SPARC: 

-  Developed  by  University  of  Karlsruhe,  Germany 

-  A  multi-block  structured  code 

-  Multigrid  strategy  is  available 

-  Parallelized  code 

-  Turbulence  models:  Spalart-Allmaras,  Baldwin-Lomax, 

Lauder  and  Sharma’s  k-r  model  and  Spezial’s  non-linear 
k-r  model 

•  LU-SSOR  implicit  scheme 

•  Artificial  dissipation  term 
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Results  and  Discussion 


•  Numerical  errors 

•  Effect  of  the  Azimuthal  Location  of  the  Flow 
Effectors 

•  Effect  of  the  Axial  Location  of  the  Flow  Effectors 

•  Effect  of  the  Height  of  the  Flow  Effectors 

•  Effect  of  the  Number  of  Flow  Effectors 

•  Comparison  with  wind  tunnel  test  data 

•  Efficiency  of  the  flow  effectors 
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Results  and  Discussion 


Effect  of  grid  distribution  for  the  body-alone  configuration 


Uniform  Mesh 

Non-uniform  Mesh 

AOA  =  15° 

CY=  0.00358 

C/Cz=  0.00280 

CY=  0.02033 

C/Cz  =  0.01288 

AO  A  =  17° 

CY=  0.00567 

Cy/Cz=  0.00335 

CY=  0.03043 

C/Cz  =  0.01455 
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Results  and  Discussion 


•  Effect  of  the  Azimuthal  Location  of  the  Flow  Effectors 
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Results  and  Discussion 


•  Effect  of  the  Axial  Location  of  the  Flow  Effectors 


Front  location: 
x/D=0.844 

Back  location: 
x/D=1 .286 


Single  FE,  AOA=17° 


Results  and  Discussion 


Effect  of  the  Height  of  the  Flow  Effectors 


FE  Height 
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1H 

1.5H 

2H 

c/cz 
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Results  and  Discussion 


•  Effect  of  the  Number  of  Flow  Effectors 


♦  RFE1-270  (spare)  o  RFE1-270  (WT) 
a  RFE3-270  (spare)  a  RFE3-270  (WT) 

•  RFE1-90  (spare)  o  RFE1-90  (WT) 

■  RFE3-90  (spare)  □  RFE3-90  (WT) 


♦  RFE1-225  (spare)  o  RFE1-225  (WT) 

▲  RFE3-225  (spare)  a  RFE3-225  (WT) 
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Results  and  Discussion 


•  Comparison  with  wind  tunnel  test  data 

-  Normal  and  axial  force  coefficients 
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Results  and  Discussion 


•  Comparison  with  wind  tunnei  test  data 

-  side  force  coefficients 
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o  RFE1-225  (WT) 
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Results  and  Discussion 


•  Comparison  with  wind  tunnei  test  data 


-  moment  coefficients 


♦  RFE1-225  (spare)  o  RFE1-225  (WT) 


AOA 


Results  and  Discussion 


•  Efficiency  of  the  fiow  effector 

-  side/normal  force  coefficients 


■  RFE1-135  (spare)  □  RFE1-135  (WT) 
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Lauder  and  Sharma’s  k-r  model 
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Conclusions 


•  The  maximum  peak  side  forces  occurred  at  AOA  =17°  for  all  the 
configurations  with  flow  effectors  at  the  front  axial  location 

x/D  =  0.844. 

•  Among  the  configurations  studied,  the  single  FE  situated  at  an 
azimuthal  angle  of  225°  was  the  most  efficient  in  generating  the  side 
forces. 

•  The  side  force  gradually  decreased  as  the  FE  position  moved  to  the 
horizontal  symmetry  plane  of  the  body. 

•  The  front  configuration  was  more  efficient  than  the  back  configuration 
in  influencing  the  vortex  system  behaviour. 

•  The  side  force  did  increase  with  the  number  of  FEs;  however,  the 
advantage  was  not  dramatic. 
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Conclusions 


•  the  side  force  coefficients  on  the  starboard  side  of  the  body  agreed 
well  with  the  wind  tunnel  data,  but  those  on  the  port  side  of  the  body, 
which  had  identical  numerical  counterparts,  were  consistently  lower 
than  the  wind  tunnel  data.  The  numerical  normal  and  axial  forces 
agreed  with  the  wind  tunnel  data  well. 

•  The  numerical  study  showed  that  the  lateral  forces  generated  by  case 
RFE1-225  were  on  average  about  25%  of  the  normal  forces.  With 
such  a  small  device  installed  at  the  nose  of  the  missile,  the 
lateral  force  gain  was  quite  remarkable.  This  re-confirms  that  the 
concept  of  using  small  devices  at  the  nose  of  the  missile  to  control  the 
lateral  forces  is  feasible. 
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4.3  Impact  of  Nose-mounted  Micro-structures  on  the  Aerodynamics  of  a 
Generic  Finned  Missile 

presented  by  D.  Corriveau,  DRDC  -  Valcartier 

Balance  force  measurements  were  performed  on  a  generic  missile  configuration  with 
finned  mounted  in  cruciform  configuration.  The  baseline  experimental  model  used  consists  in  a 
finned  missile  having  a  conical  nose  and  an  aspect  ratio  (L/D)  of  13.  Tests  were  performed  in  the 
DRDC  indraft  wind  tunnel  at  a  free  stream  Mach  number  of  1.5  and  for  angles-of-attack  varying 
between  0.0°  and  20.0°.  Several  configurations  of  the  nose-mounted  flow  effectors  were 
investigated.  The  results  show  that  no  peak  in  the  side  force  generated  is  attained  for  angles  of 
attack  up  to  20.0°.  It  was  also  found  that  the  side  force  is  higher  when  the  flow  effectors  are 
located  closer  to  the  nose  tip.  By  indexing  the  flow  effectors  to  various  roll  positions  on  the  nose, 
it  was  determined  that  the  maximum  side  force  occurs  at  an  angular  position  of  225°.  On  the 
downside,  for  missile  equipped  with  fins,  the  side  force  generated  by  the  flow  effectors  is 
relatively  small  in  comparison  to  the  normal  force. 
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Objectives 

•  Determine:  Optimum  flow  effector  configurations 

•  Determine:  Range  of  AOA  for  which  a  side  force 

is  generated 

•  Understand:  Mechanisms  through  which  the 

flow  effectors  impact  on  the  side 
force’s  magnitude 

•  Evaluate:  Use  of  vortex-induced  side  force  to 

achieve  yaw  control 
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ruy  Background 

Slender  Body  Vortex  Shedding 


!>:*»- ill 


•  Vortices  form  on  the  missile’s  leeside 
operating  at  AOA 

•  Early  studies  dates  back  to  the  50’s 

•  Early  concerns  related  to  vortex 
interactions  with  the  control  surfaces 
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Background 

Side  Force  Control  Using  Micro-Structures 


AFRL  MN  (2003) 

>  Air  Force  Basic  Research  Shape 

>  A/B  range  tests 

>  Wind  tunnel  test  DRDC  Valcartier 


Orbital  Research  /  DARPA  (2002) 

>  High  AOA  (+30°  to +60°) 

>  Low  aspect  ratio  model  (I/d  =  4.5) 

>  Incompressible  subsonic  flow 

>  Control  existing  side  force 


Leu  et  al.  (2005) 

>  Low  AOA  (0°  to +60°) 

>  Low  aspect  ratio  model  (I/d  =  5.0) 

>  Incompressible  subsonic  flow 

>  Control  existing  side  force 
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f  \$J  Background 

Side  Force  Control  Using  Micro-Structures 


Current  Study  (2004) 

>  Low  AOA  (0°  to  +20°) 

>  Medium  aspect  ratio  model  (I/d  =  13) 

>  Supersonic  flow  (M  =  1 .5) 

>  Generate  side  force  for  yaw  control 


Defence  R&D  Canada  -Valcartier  •  R  &  D  pour  la  defense  Canada  -  Valcartier 


rto  Presentation  Overview 


Objectives 

Background 

Test  setup  and  instrumentation 
Missile  configuration 
Experimental  results 
Conclusions 


Defence  R&D  Canada  -  Valcartier  •  R&D  pour  la  defense  Canada  -  Valcartier 


4.3-5 


r^y  Test  Setup  and  Instrumentation 


DRDC  Trisonic  Wind  Tunnel 


950  m3 

Air  Storage  Bag 


950  mJ 

Air  Storage  Bag 


•  Indraft  wind  tunnel 

•  Test  section:  2’x  2’ 

•  0.2  <  Ma  <  4.0 

•  -20°  <  AOA  <  +20  ° 

Defence  R&D  Canada  -  Valcartier 


•  Test  duration:  5-11  sec 

•  Turnaround  time:  30  min. 
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r^y  Missile  Configuration 

Baseline  Geometry  With  Conical  Nose 


Aspect  Ratio  L/D  =  13.0 
Nose  Aspect  Ratio  Ln/D  =  3.0 
Conical  Nose 

4  Fins  in  +  Configurations  (Removable) 
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Missile  Configuration 

Baseline  Geometry  With  Conical  Nose 


Missile  Configuration 

Baseline  Geometry  With  Conical  Nose 
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f  ^  Missile  Configuration 

Nose  Arrangement 
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Missile  Configuration 

Flow  Effector  Geometries 
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f  \i_)  Missile  Configuration 

Flow  Effector  Configurations 
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Experimental  Results 

Lateral  Force  Coefficient  vs  AOA 


•  4  fins  in  +  configuration 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  270 


•  Front  row 

•  Ma=  1.5 
•Re/m=1 5.2x1 06 
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Experimental  Results 

Lateral  Force  Coefficient  vs  AOA 


Defence  R&D  Canada  -  Valcartier 


•  4  fins  in  +  configuration 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  225° 

•  Front  row 

•  Ma=  1.5 
•Re/m=1 5.2x1 06 
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r^y  Experimental  Results 

Lateral  Force  Coefficient  vs  AOA 


•  4  fins  in  +  configurations 

•  Conical  nose  LN/D=  3.0 

•  Flow  effector:  1  key 

•  Front  row 

•Angular  position:  0°  -  360° 

•  Ma=  1.5 

•  Re/m=  1 5.2x1 06 
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Experimental  Results 

Lateral  Force  to  Normal  Force  Coefficient  Ratio  vs  AOA 
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•  4  fins  in  +  configuration 

•  Conical  nose  LN/D=  3.0 

•  Flow  effector:  1  key 

•  Angular  position:  0°  -  360° 

•  Front  row 

•  Ma=  1.5 

•  Re/m=  15.2x1 06 
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Lateral  Force  Coefficient  vs  AOA 
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Experimental  Results 

Axial  Force  Coefficient  vs  AOA 
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•  4  fins  in  +  configuration 

•  Conical  nose  LN/D=  3.0 
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•  Front  row 
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Axial  Force  Coefficient  vs  AOA 


•  4  fins  in  +  configuration 

•  Conical  nose  LN/D=  3.0 

•  Flow  effectors  centered  at  270° 

•  Front  row 
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Normal  Force  Coefficient  vs  AOA 
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Normal  Force  Coefficient  vs  AOA 
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•  Flow  effectors  centered  at  270 
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f  Conclusions 


Peak  Side  force  magnitude  measured  between  AOA  of 
15.0°  and  20.0° 

Locating  the  flow  effectors  closer  to  the  nose  tip  results 
in  higher  side  forces 

Maximum  side  forces  measured  at  an  angular  position 
of  225° 

Side  forces  generated  by  the  micro-structures  are 
symmetric  for  complementary  flow  effectors  angular 
position 

Apart  from  side  force,  the  nose-mounted  micro¬ 
structures  have  little  impacts  on  the  other  aerodynamic 
coefficients 
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4.4  Mini-scale  PIV  Measurements  during  Unsteady  Excitation  of  Mini- 
Leading  Edge  Flaps  Used  for  Separation  Control  and  Development  of  a  Bi¬ 
fold  Five  Component  Force  Balance 

presented  by  X-Z.  Huang,  IAR/NRC 

Reduction  of  boundary  separation  is  a  key  issue  to  increase  lift  and  control  power 
especially  at  take-off,  landing  and  post-stall  flight  conditions.  Leading-edge  flaps,  either  static  or 
dynamic,  are  one  of  the  promising  methods  that  has  been  studied.  Miniature  flaps  activated  by 
micro-actuators  may  reduce  power  demands  of  such  devices.  A  clear  understanding  on  how 
aerodynamic  forces  can  be  altered  by  manipulating  the  micro-flow  effectors  is  a  prerequisite  for 
its  application.  To  this  end,  an  unsteady  excitation  of  a  micro-leading-edge  flap  concept  was 
introduced  and  studied  at  IAR.  The  experiments  were  conducted  in  the  IAR  1520  water  tunnel 
with  a  two  dimensional  wing  profile  (NACA  0012)  as  test  model.  Mini-scale  PIV  measurements 
showed  that  when  the  leading-edge  was  oscillating,  the  boundary  layer  separation  area  was 
reduced  dramatically.  Parameter  studies  found  that  when  the  non-dimensional  frequency,  fc/  U  oo 
>1,  the  effectiveness  becomes  smaller  as  the  frequency  increased.  An  analysis  shows  that  the 
effectiveness  of  the  moving  surface  is  related  to  the  increased  viscosity  present  in  the  sub-layer 
with  the  increment  in  effectiveness  being  related  to  the  frequency  of  the  moving  wall. 

In  a  connected  activity,  a  five-component,  bi-fold  force  balance  was  developed  for  data 
measurement  in  water/wind  tunnel  open  test  sections.  Since  the  two  dimensional  profile  wing  is  a 
half-model  where  only  the  wing  is  immersed  in  the  water  tunnel,  a  short  balance  is  needed 
because  the  center  pressure  of  the  model  should  be  as  close  to  the  pivot  of  the  motion  control 
system  as  possible.  On  the  other  hand,  requirements  for  sensitivity,  linearity  and  gauge  space 
does  not  allow  the  balance  elements  to  be  too  short.  The  bi-fold  concept  keeps  the  length  of  each 
element  as  long  as  possible  while  shortening  the  total  length  of  the  balance.  Numerical 
calculations  based  on  the  finite  element  method  were  conducted  and  compared  with  data 
measured  from  a  prototype  bi-fold  balance.  The  experimental  results  confirmed  that  the  concept 
could  be  successfully  used  in  half-model  experiments. 
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Topic  1 

Development  of  a  New  Concept  of  Bi-Fold , 
Five-Component,  Half-Model  Balance  used  in 
IAR  Water  Tunnel 

♦>  Introduction 

❖  Principle  of  bi-fold  balance 

❖  Geometrical  description  of  initial  balance 

❖  Gauges  and  bridges 

❖  Calibration  results 

❖  Further  configuration  studies  and 
recommendations 
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bi-fold  balance 


Forward  elements 


Rolling  elements 


aft  elements 


Bi-fold  balance 


Normal  balance 
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Attach  to  model 


4.4-2 


NRCa  eraspa  Ec.cout 


4.4-3 


NRCa  eraspa  Ec.cout 


Further  configuration 
studies  and 
recommendations 


.150  — 
0.100— ^ 

3$01_ 

TZ 

1 

.220  J 

— 

R 

K 

.080 

.150 


4.4-4 


Topic  2 

Development  of  Mini-Scale  PIV 
Methodology  and  its  Applications 
on  Wall-Bounded  Flows 


♦>  Introduction 


♦>  Experimental  set-up 
♦>  Experimental  results 
♦>  Analysis  and  conclusions 


Development  of  Mini-Scale  PIV 
Methodology  and  its  Applications 
on  Wall-Bounded  Flows 


Introduction 

Wall-bounded  flow  is  full  of  vorticities 

Aerodynamic  behavior  is  mainly  determined  by  vorticity 
generation  in  boundary  layer  and  convection  downstream 

Leading-edge  area  is  crucial  for  voticity  generation  and 
convection  downstream 

The  effectiveness  of  mini  flow  effectors  depends  on  the  studying 
of  the  mini-scale  wall-bounded  flow,  especially  at  leading-edge 
area 


4.4-5 


ulence  structure  at  various  heights 
in  a  low-Re  boundary  layer  condition 
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Top  and  side  views  of  wall-bounded  flow 
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I ini-scale  PIV  measurements  on  wall- 
bounded  flow  at  the  mid  chord  of 
NACA  0012,  Uao=2.7in/sec,  Re=13450 


rini-scale  PIV  measurements  of  wall- 
bounded  flow  on  the  mid  chord  of 
NACA  0012,  U°o=5.3in/sec,  Re=26900 
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ini-scale  PIV  measurements  of  wall- 
bounded  flow  on  the  mid  chord  of 
NACA  0012,  Uoa=7.9in/sec,  Re=40350 
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lini-scale  PIV  measurements  of  wall- 
bounded  flow  at  the  leading-edge  of 
NACA  0018,  U°o=2. 7 in/sec,  Re=13450 
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♦>  The  resolution  could  reach  50|um  and  view  area  could  be  2~5  mm 


❖  The  instability  in  the  open  shear  flow  may  be  classified  as  either  noise 
amplifiers  with  extrinsic  dynamics  or  as  oscillators  with  intrinsic 
dynamics 

❖  For  the  wall-bounded  flows,  the  velocity  gradient  at  the  wall  is  the 
highest  one  resulting  in  the  maximum  production  of  turbulent  energy 

❖  That  energy  source  is  exported  to  the  rest  of  the  flow  as  a  spatial  energy 
flux.  When  the  flow  velocity  is  very  low  and  the  outside  flow  is  in 
laminar  state,  the  outside  flow  acts  as  energy  sinks.  The  dissipation 
process  will  absorb  the  energy  flux.  At  this  situation  the  flow  acts  as  a 
noise  amplifiers  and  the  flow  will  be  back  to  laminar  eventually  as 
observed  in  the  experiments. 

❖  If  the  flow  velocity  is  high  enough  and  the  flow  state  crosses  the  border 
of  the  oscillators  with  intrinsic  dynamics ,  those  small  eddies  may 
trigger  large  flow  instability 


Boundary  layer  separation  control 
by  unsteady  excitation  of  mini-scale 
leading-edge  flaps 
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Topic  3 

oundary  layer  separation  control 
by  unsteady  excitation  of  mini- 
scale  leading-edge  flaps 


Smoke  flow  visualization  on 
different  leading-edge  flap  angles 
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Mini-scale  leading-edge  flap 


boundary  layer  separation  control 
by  unsteady  excitation  of  mini¬ 
scale  leading-edge  flaps 
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oundary  layer  separation  control 
by  unsteady  excitation  of  mini¬ 
scale  leading-edge  flaps 


Mini-scale  leading-edge  flap  inside  (static) 


oundary  layer  separation  control 
by  unsteady  excitation  of  mini¬ 
scale  leading-edge  flaps 


Mini-scale  leading-edge  flap  outside  (static) 
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oundary  layer  separation  control 
by  unsteady  excitation  of  mini¬ 
scale  leading-edge  flaps 


Mini-scale  leading-edge  flap  oscillating  at  0.5  Hz 


oundary  layer  separation  control 
by  unsteady  excitation  of  mini¬ 
scale  leading-edge  flaps 


Mini-scale  leading-edge  flap  oscillating  at  1  Hz 
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rini-scale  PIV  measurement  for  mini¬ 
scale  leading-edge  flap  inside  the  wing 

profile  (static) 
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ini-scale  PIV  measurement  for  mini¬ 
scale  leading-edge  flap  outside  the  wing 

profile  (static) 
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Future  work 
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Conclusion 


1 .  New  concept  of  bi-fold  balance  has  been 
developed  and  constructed 

2.  Mini-scale  PIV  method  can  be  used  to 
visualize  wall-bounded  flow  field 

3.  Unsteady  excitation  of  mini-scale  leading- 
edge  flaps  are  able  to  create  enough  control 
force  on  boundary  layer  separation 

4.  Full  model  experiments  remain  to  be 
completed 
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5.  Microactuator  Development 


5.1  SMA  Antagonistic  Actuator  Modeling  and  Validation  Results 

presented  by  O.  Boissonneault,  Numerica 

In  the  field  of  micro-machines,  shape  memory  alloys  (SMA)  are  receiving  increased 
attention  because  of  their  high  power-to-volume  ratio  and  large  recoverable  deformation. 
However,  these  materials  have  a  complex  thermomechanical  behavior  due  to  the  hysteresis 
related  to  stress,  strain  and  temperature.  A  constitutive  model  that  can  consider  thermal  or 
mechanical  loads  or  a  combination  of  loads  was  developed  in  order  to  predict  shape  memory 
effect  as  well  as  superelasticity. 

The  parameters  for  the  constitutive  model  in  the  form  of  a  phase  diagram  is  required. 
They  were  measured  using  a  0.1  mm  diameter  SMA  wire  with  a  Rheometrics  (RSA-II)  Dynamic 
Mechanical  Analyser  in  a  series  of  isothermal  tensile  tests.  All  tests  were  in  the  range  of  30°C  to 
130°C. 

A  small  test  bench  was  built  to  compare  the  predicted  and  experimental  results.  SMA 
wire  measuring  approximately  100  mm  in  length  was  excited  by  short  and  intense  voltage  pulses 
measuring  between  4.5  to  10  volts.  The  contraction  of  the  SMA  was  used  to  lift  a  constant  150  g 
load. 

An  another  test  bench  was  built  to  study  the  performance  of  SMA  wire  in  an  antagonistic 
set-up.  A  series  of  2  or  3  volt  step  inputs  were  alternately  injected  into  each  wire  to  characterize 
the  peak-to-peak  displacement  and  the  motion  time  constant.  An  antagonistic  actuator  model 
based  on  the  hybrid  SMA  model  predicted  reasonably  well  the  displacement-time  results. 
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SMA  Actuator  Modeling  and  Validation 

Results 

Olivier  Boissonneault,  RDDC-NUMERICA 


Action  list 


8  -  Characterize  Dynalloy  SMA  wire  and  construct  a 

small-scale  test  bench 

9  -  Examine  ability  of  SMA  model  to  predict 

response  under  complex  force  profile 

10  -  Investigate  methods  to  increase  frequency 

response  of  the  SMA  wire  given  force 
requirements 

Develop  method  to  analyse  compliant  beam. 
Integrate  method  with  SMA  model  to  determine 
appropriate  size  for  beam  elements. 

Integrate  control  law  in  test  bench  environment. 
Explore  ability  to  track  command  signal. 
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Characterization  of  the  SMA  wire 
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Results  for  4.5  volts 
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Results  for  10  volts 
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,v 


Displacement  of  1  mm 


Voltage 

(V) 

Heating 

(ms) 

Cooling 

(ms) 

Total  Time 

(ms) 

Oscillations 

Power 

(W) 

Energy 

(J) 

10 

26 

1070 

1096 

+  +  +  + 

10 

0.26 

8 

40 

1210 

1250 

+  +  + 

6.4 

0.256 

6 

100 

1160 

1260 

+  + 

3.6 

0.36 

4.5 

175 

1140 

1315 

+ 

2.03 

0.304 

P  =  E7 
R 

*  Hypothesis  :  the  electrical  resistance  is  10  ohms 


Action  list 


8  -  Characterize  Dynalloy  SMA  wire  and  construct  a 

small-scale  test  bench 

9  -  Examine  ability  of  SMA  model  to  predict 

response  under  complex  force  profile 

10  -  Investigate  methods  to  increase  frequency 

response  of  the  SMA  wire  given  force 
requirements 

13  -  Develop  method  to  analyse  compliant  beam. 
Integrate  method  with  SMA  model  to  determine 
appropriate  size  for  beam  elements. 

Integrate  control  law  in  test  bench  environment. 
Explore  ability  to  track  command  signal. 
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Antagonistic  test  bench 
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Modification  of  the  model  hybrid 
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Action  list 


8  -  Characterize  Dynalloy  SMA  wire  and  construct  a 

small-scale  test  bench 

9  -  Examine  ability  of  SMA  model  to  predict 

response  under  complex  force  profile 

10  -  Investigate  methods  to  increase  frequency 

response  of  the  SMA  wire  given  force 
requirements 

Develop  method  to  analyse  compliant  beam. 
Integrate  method  with  SMA  model  to  determine 
appropriate  size  for  beam  elements. 

Integrate  control  law  in  test  bench  environment. 
Explore  ability  to  track  command  signal. 
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Future  works 


•  Calibration  of  flow  effectors  strain  gage  (last 
generation) 

•  System  identification  of  the  actuators  for  control 

•  Validation  of  the  closed  loop  control 
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5.2  Compliant  Flow  Effector  Design  for  Missile  Side  Force  Control 

presented  by  F.  Wong,  DRDC  -  Valcartier 

Active  flow  control  can  be  used  in  a  control  actuation  system  to  generate  control 
authority  for  supersonic  missiles.  Placing  small  flow  effectors  close  to  the  nose  tip  permits 
modulation  of  vortex  shedding  to  generate  controllable  side  forces  and  moments.  The  lack  of 
volume  for  an  actuator  in  the  missle  nose  necessitates  a  different  approach  be  taken  for  its  design. 
To  this  end,  an  actuator  comprised  of  a  shape  memory  alloy  actuator  (SMA)  coupled  to  a 
compliant  mechanism  was  conceived.  The  mechanism  acts  as  a  transmission  to  convert  the  force- 
displacment  characteristics  of  the  SMA  wire  to  the  displacement  characteristics  required  to 
generate  control  forces. 

Three  compliant  mechanisms  were  studied.  Pivot,  Trapezoid  and  Inverted  Trapezoid 
geometries  were  selected  for  their  simplicity.  The  main  constraints  for  a  viable  mechanism 
geometry  were  maximizing  the  tip  displacement  while  keeping  the  post  base  stresses  below  the 
structure’s  yield  stress.  It  was  found  that  a  trapezoid  mechanism  gave  a  tip  displacement  of  870 
micron  while  remaining  within  the  yield  stress  limits  of  the  material.  The  Pivot  and  Inverted 
Trapezoid  geometries  could  not  transform  the  SMA  force  and  displacement  input  without 
yielding  the  structure. 
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f  \i_)  Background 


Wind  tunnel  studies  on  finless  missile  showed 
that  the  magnitude  of  the  side  force  can  be 
modulated  by  positioning  micro-flow 
effectors  at  different  angular  positions  around 
central  axis. 

Shape  memory  alloy  study  developed  a 
constitutive  model  that  captures  the  temporal 
thermomechanical  behavior  of  the  material. 

Control  algorithm  was  synthesized  to 
compensate  for  hysteretic  behaviour  and 
provide  good  tracking  performance  at 
moderately  high  frequencies. 
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Compliant  Mechanism  Requirements 


Defence  R&D  Canada  -  Valcartier  •  R  &  D  pour  la  defense  Canada  -  Valcartier 


Flow  effectors:  1  mm  x  2.67  mm  x  0.76  mm 
Diameter  at  1st  row:  10  mm 
Operating  parameters: 

a)  1  mm  flow  effector  deflection 

b)  150-200  g  SMA  applied  force 

c)  no  yielding  of  material 

d)  1  Hz  operation 


Problem  Statement 


Seek  a  structure  that  will 
meet  kinematic  requirements 
while  at  the  same  time  meet 
mechanical  requirements. 

Too  flexible  a  structure  will 
not  support  loads. 

Too  rigid  a  structure  will 
require  high  forces  to  achieve 
motion. 

Synthesis  procedure  finds  an 
optimum  structure  that 
satisfies  both  requirements. 
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r\,y  Synthesis  Criteria 


Required  kinematic  motion 
Available  design  space 
Desired  geometric  or  mechanical  advantage 
Material  properties 
Stress  or  strain  limitations 


Current  synthesis 
methods  address 
these  criteria. 


Buckling  instabilities 
Weight  limitations 
Dynamic  behaviour 
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3505.0 

L2LEN 

(micron) 

8839.0 

L3HGT 

(micron) 

2540.0 

L4HGT 

(micron) 

1499.0 

L5HGT 

(micron) 

2032.0 

ARMWID 

(micron) 

508.0 

LEVHGT 

(micron) 

508.0 

TOPWID 

(micron) 

203.0 

BASWID 

(micron) 

1575.0 

FEWID 

(micron) 

2794  0 

LEVLEN 

(micron) 

12344.0 

FETHK 

(micron) 

787.0 

THGT 

(micron) 

1016.0 

PIVRAD 

(micron) 

191.0 

SY  allow 

(MPa) 

497.0 

Modulus 

(MPa) 

72.0 

FSMA 

isH 

150.0 
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>v 


f  OutputY 


1 

2 

3 

TOPWID 

508.0 

(micron) 

787.4 

(micron) 

787.4 

(micron) 

BASWID 

787.4 

(micron) 

508.0 

(micron) 

254.0 

(micron) 

InputX 

21.3 

(micron) 

22.9 

(micron) 

56.4 

(micron) 

OutputY 

41.3 

(micron) 

44.3 

(micron) 

154.8 

(micron) 

BaseSY 

57.5 

(MPa) 

167.8 

(MPa) 

713.7 

(MPa) 

TipSY 

54.9 

(MPa) 

16.1 

(MPa) 

14.1 

(MPa) 

Out/In  Disp 

1.9 

1.9 

2.7 

BaseSY/SY 

0.12 

0.34 

1.44 

TipSY/SY 

0.11 

0.03 

0.03 

L1LEN 

L2LEN 

L3HGT 

L4HGT 

L5HGT 


TOPWID 

BASWID 

ARMWID 

LEVHGT 

FEWID 

LEVLEN 

FETHK 

THGT 


SY  allow 
Modulus 
FSMA 
FAERO 


(MPa) 

(MPa) 


(micron) 

(micron) 

(micron) 

(micron) 

(micron) 


3759.0 

8839.0 

2540.0 

1498.6 

2032.0 


(micron) 

(micron) 

(micron)^ 

(micron) 

(micron) 

(micron) 

(micron) 

(micron) 


508.0 

787.0 

762.0 

508.0 

2794.0 

14760.0 

787.0 

1016.0 


497.0 

150.0 

150.0 

0.172 
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Ra7  Trapezoidal  Post 
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LI  LEN 

(micron) 

3759.0 

L2LEN 

(micron) 

8839.0 

L3HGT 

(micron) 

2540.0 

L4HGT 

(micron) 

1498.6 

L5HGT 

(micron) 

2032.0 

TOPWID 

(micron) 

508.0 

BASWID 

(micron)  I 

787.0 

ARMWID 

(micron) 

762.0 

LEVHGT 

(micron) 

508.0 

FEWID 

(micron) 

2794.0 

LEVLEN 

(micron) 

14760.0 

FETHK 

(micron) 

787.0 

THGT 

(micron) 

1016.0 

SY  allow 

(MPa) 

497.0 

Modulus 

(MPa) 

150.0 

FSMA 

M _ 

150.0 

FAERO 

W 

0.172 
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5.2-4 


Validation  of  Analysis 


Measured  flow  effector  tip  and  T-section 
displacements  and  compliant  link  strains. 
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Displacement  Results 
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rxo  Compliant  Link  Strain  Results 


Force  (g) 
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.V 


Final  Configuration 
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RD7  Final  Flow  Effector 


L1LEN 

L2LEN 

L3HGT 

L4HGT 

L5HGT 

(micron)  3607.0 

(micron)  21158.0 

(micron)  3429.0 

(micron)  2540.0 

(micron)  2286.0 

TOPWID 

BASWID 

ARMWID 

LEVHGT 

FEWID 

LEVLEN 

FETHK 

THGT 

(micron)  254.0 

(micron)  279.0 

(micron)  762.0 

(micron)  762.0 

(micron)  2667.0 

(micron)  24765.0 

(micron)  787.0 

(micron)  1207.0 

SY  allow 

(MPa)  882.0 

reometry 


FEWID  -I  TOPWID 


1 


-|  BASWID  |- 


FSMATop 

(9) 

150.0 

0.0 

150.0 

0.0 

185.0 

185.0 

(9) 

FSMA  Bot 

(9) 

0.0 

150.0 

0.0 

150.0 

0.0 

150.0 

(9) 

FAero 

(N) 

0.172 

0.172 

0.0 

0.0 

0.172 

0.172 

(N) 

SMA  Top  X 

(micron) 

154.7 

-189.2 

367.2 

23.3 

240.4 

263.6 

(micron) 

SMA  Top  Y 

(micron) 

-136.2 

202.7 

-345.1 

-6.2 

-216.7 

-222.9 

(micron) 

SMA  Bot  X 

(micron) 

16.7 

14.5 

23.3 

21.1 

22.2 

43.2 

(micron) 

SMA  Bot  Y 

(micron) 

-135.9 

203.0 

-344.8 

-5.9 

-216.3 

-222.3 

(micron) 

FE  Frt  X 

(micron) 

51.8 

-40.9 

112.1 

19.4 

78.0 

97.4 

(micron) 

FE  Frt  Y 

(micron) 

662.4 

-1351.6 

2078.1 

64.2 

1147.3 

1211.4 

(micron) 

FE  Bek  X 

(micron) 

60.7 

-61.4 

142.4 

20.3 

93.9 

114.2 

(micron) 

FE  Bek  Y 

(micron) 

586.5 

-1175.1 

1817.7 

56.1 

1010.6 

1066.7 

(micron) 

Link  Tip  X 

(micron) 

59.7 

-49.2 

128.8 

19.9 

89.8 

109.6 

(micron) 

Link  Tip  Y 

(micron) 

4.7 

-6.8 

12.1 

0.5 

7.6 

8.1 

(micron) 

Link  Base  SY 

(MPa) 

449.0 

-100.5 

798.8 

249.3 

635.4 

884.6 

(MPa) 

Link  Tip  SY 

(MPa) 

19.1 

-481.8 

344.4 

-156.4 

99.4 

-57.0 

(MPa) 

Link  Base  ey 

(--) 

0.0032 

-0.0012 

0.0059 

0.0016 

0.0045 

0.0061 

(--) 

Out/In  Disp 

4.3 

7.1 

5.7 

2.8 

4.8 

4.6 
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Summary 


Compliant  mechanism  principles  used  to  develop  a 
transmission  mechanism  to  transform  the  force- 
displacement  characteristics  of  the  actuator  to  the 
required  force-displacement  characteristics  of  the 
flow  effector. 

Initial  mechanism  produced  90%  of  the  desired 
displacement  without  yielding  the  material. 

Design  validated  with  displacement  and  strain  data 
under  applied  forces. 

Design  revised  to  produce  1  mm  displacement  and 
adjusted  for  missile  model  envelope  constraints. 

Design  of  four  flow  effector  assembly  produced. 
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6.  Micro-Fabrication  Technology 


6.1  Precision  Fabrication  of  Active  Flow  Effectors  for  Supersonic  Missile 
Control 

presented  by  S.K.  Nikumb,  IMTI/NRC 

Development  of  highly  precise  flow  effectors  based  on  the  results  of  studies  in 
aerodynamics,  control  system  requirements  and  design  specifications,  posed  significant 
challenges  to  manufacture  them.  In  this  presentation,  we  report  our  results  on  the  precision 
fabrication  of  the  test  bench  and  final  optimum  geometries  of  the  flow  effectors.  To  achieve 
desired  geometric  tolerances,  repeatability  and  good  surface  finish,  different  machining  methods 
were  employed  and  their  relative  performance  was  evaluated.  The  machining  method  was 
optimized  during  the  fabrication  of  four  versions  of  the  flow  effectors  made  from  Aluminum, 
Brass,  Shape  Memory  Alloy  (SMA)  and  Titanium  materials.  Laser  machining,  micro-hole 
drilling,  laser-milling  as  hybrid  process  and  micro  milling  processes  were  used  in  the  fabrication 
of  the  flow  effectors.  Based  on  the  results  obtained  from  the  initial  bench  tests,  several  samples  of 
Titanium  flow  effectors  were  successfully  fabricated  with  high  dimensional  accuracy,  <6  micron 
tolerance  and  with  a  cut  edge  surface  finish  down  to  several  tens  of  nanometers.  The  samples 
were  measured  and  characterized  using  optical  inspection  techniques  SEM,  Microscopes,  WYKO 
and  shadow  graph  instruments.  Interesting  observations  were  made  on  tool  wear  and  its  effect  on 
the  cut  edge  surface  finish. 
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Precision  Fabrication  of  Flow  Effectors 
for  Supersonic  Missile  Control 


SUWAS  NIKUMB 
NRC-IMTI 

Suwas.  nikumb@nrc.  ca 


Ij.!  National  Research  Conseil  national  O  Y^l  SI  Z"’]  Si 

■  Council  Canada  de  recherches  Canada  X^CXJL  JLCLVACX. 

DRDC  Final  Meeting  April  19th,  2006 


Introduction 


Precision  Fabrication  of  the  Flow  Effector 

Objective:  To  achieve  high  dimensional  tolerances,  repeatability,  and 
good  edge  surface  finish 

Version  1:  Process  selection,  Parameters  development,  Laser  cutting 
of  flow  effectors  out  of  Al,  SMA,  Ti,  and  Brass,  cut  and  drill 
quality 

Version  2:  New  geometry,  Experimental  Verification,  Fabrication  of 
flow  effectors  out  of  Al,  SMA,  Ti,  and  Brass  on  materials 
Version  3:  Dimensional  changes,  Test  bed  set  up 
Version  4:  Final  version  fabrication  in  Titanium,  Challenges,  WYKO 
Measurements,  Fixture  development,  Surface  finish,  Tool 
wear,  Observations  with  and  without  coolant. 


6.1-1 


-0--E3- 


1  “TC 

p 

f 

RAD  pour  li  defense  Canada  -  Vak*rtoer 

QLGBEC,  CttAEA  GSJtfS 

_ 

0 

secti^Vrmes 

FLIGHT  OVNNMCS  GROUP 

PRECIS  ON  WEAPONS  SECTION 

M1E 

ll| 

SI  f 

CUMr’LlAN  1  btAM  UlbX  VLKbKJN  Itbl  A 1  LLItK 

^05042801 DFT 

FW0504280I 

SCAIE 

I/I 

OE 

B 

Geometric  Specs  determination,  Post  thickness,  250  (am  supports  requirements,  processes  -  Laser,  Micro-milling  and  Both,  Samples  prepared 
and  evaluated.  Laser  Process  too  long.  Materials  Brass,  Aluminum  7075,  A  number  of  trials  conducted  however  samples  were  never  tested. 


Al  7075  Edge  and  drill 


x  80 


x  32 


Micro  milling  /  laser 
combined  process 
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Micro  milling  /  laser 
combined  process 


Flow  Effector  Version  2 
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FLOW  EFFECTOR  V2 

FW05042808DFT  FW05042808 


New  geometry  with  triangular  post,  Samples  for  all  materials  prepared  and  tested/no  consistency  in  dimensions  of  diff  materials,  Holes  were 
made  but  laser  drill  was  removing  post  material,  conical  entrance  small  exit  debris  not  removed  fully,  shielding  required  below. 
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Flow  Effector  Version  3 
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Conceptual  version  with  several  dimensional  changes,  measuring  improvements.  Never  made  samples  and  moved  on  to  version  4 
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Flow  Effector  Version  4 
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Main  change  was  in  length  X2,  angular  shape,  difficulty  in  measuring  tolerances.  Materials  Ti  and  SMA.  SMA  was  too  t 
Very  tight  tolerances  <0.0004”,  post  with  slight  angle  and  then  straight,  Process  used:  Micromilling 


_ Itfl CCjtftt 

integrated 

Manufacturing 

Technologies 

Institute 


Fabrication  Challenges 


Fabrication  Challenges: 

•  high  accuracy:  +0.0004”  (+6.36  (jm) 

•  high  dimensional  aspect  ratio:  0.01 1”:0. 965”  (1:87) 

•  high  dimensional  repeatability  between  parts 

•  low  surface  roughness 

•  vertically  of  machined  sides 
Micro-milling  challenges: 

•  machinability  of  Ti  /  new  “zero”  /  optimal  process  parameters  /  tool  wear 
General: 

•  Material  Hardness  -  RC  23 

•  10  +  2  parts  Micro  machined  within  tolerances 

•  Surface  finish  -  along  cutting  direction  Ra  <  600nm 

across  cutting  direction  Ra  <  800nm 

Successfully  made  10+2  samples  (in  all  19  tested) 
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Total  number  of  fabricated  Flow  Effectors:  19 


Measured  Dimensions 
Flow  Effector 


Flow  Effector  Ver.4  (F 


Dim  j 

Part# 

Req 

6 

Mill  2 . 

7 

Miii  2 . 

8 

. Mill  2 . 

9 

. Miii  2 . 

3EZ 

Mill  2 . 

11 

. Miii  2 

. 12 . 

Mill  2 . 

1 

. Miii  2 . 

23iZ 

Mill  2 

15 

. Miii  2 . 

16 

Mill  2 . 

Method 

155.00 

154.75 

154.75 

154.95 

154.85 

154.75 

154.65 

154.70 

154.90 

154.95 

154.45 

154.90 

Microscope 

i* 

. +/-i'2o 

. -0.25 

. -025 

. -0.05 

. -015 

. -025 

. -025 

. -022 

. -616’ 

. -025 

. -055 

. -010 

30.00 

30.15 

30.50 

30.50 

30.50 

30.25 

30.00 

30.25 

30.15 

30.45 

30.05 

30.15 

Microscope 

r" 

. +220 

. 0.15 

. 050 

. 050 

. 6.50 

. 025 

. 020 

. 025 

. 015 

. 045 

. 025 

. 015 

10.00 

10.20 

10.35 

10.35 

10.30 

10.25 

10.15 

10.30 

10.30 

10.35 

10.30 

10.25 

Microscope 

r 

. +040 

. 0.20 

. 025 

. 025 

. 020 

. 025 

. 015 

. 626 

. 020 

. 025 

. 020 

. 025 

„  11.00 

11.20 

11.35 

11.35 

11.35 

11.35 

11.05 

11.40 

11.20 

11.35 

11.35 

11.15 

Microscope 

+0.40 

. 020 

. 025 

. 025 

. 025 

. 025 

. 025 

. 040 

. 020 

. 025 

. 025 

. 015 

5  L 

5.00 

5.38 

5.37 

5.23 

5.14 

5.22 

5.30 

5.42 

5.13 

5.33 

5.02 

4.97 

Vision  Gauge 

Rad 

. +/-6"s 

. 0.38 

. 027 

. 023 

. 014 

. 022 

. 020 

. 042 

. 013 

. 023 

. 022 

-023 

6  R 

5.00 

5.03 

5.04 

5.06 

5.11 

5.07 

5.24 

5.15 

5.05 

5.03 

4.95 

4.96 

Vision  Gauge 

Rad 

. +/-05 

. 0.03 

. 024 

. 0.06 

. Oil" 

. 027 

. 024 

. 015 

. 025 

. 023 

. -025 

. -024 

30.00 

29.40 

29.90 

29.90 

29.60 

29.55 

29.40 

29.65 

29.60 

29.70 

29.65 

30.05 

Microscope 

...... 

+/-520 

. -0.60 

-010 

. -010 

. -040 

-045 

. -020 

-025 

. -040 

-020 

. -025 

025 

.  65.00 

64.40 

64.55 

64.50 

64.30 

64.40 

64.50 

64.45 

64.50 

64.20 

64.60 

64.60 

Microscope 

8  +/-5.60 

-0.60 

. -045 

. -050 

-070 

. -020 

-050 

. -055 

. -050 

-626 

. -040 

-040 

105.00 

106.25 

106.20 

106.50 

106.15 

106.35 

106.50 

105.60 

106.10 

106.10 

106.00 

105.70 

Microscope 

+/-5.00 

. 1 .25 

. 120 

. 150 

. 115 

. 125 

. 150 

. 020 

. 110 

. 140 

. 120 

. 070 

90.00 

92.50 

89.50 

89.60 

89.50 

89.30 

89.40 

88.95 

89.20 

89.50 

89.20 

88.90 

Microscope 

. 

+/-520 

. 250 

. -050 

. -040 

. -050 

. -070 

. -626' 

. -125 

. -020 

. -050 

. -020 

. -110 

Angle 

99.46 

99.60 

99.10 

99.50 

99.10 

99.20 

99.50 

99.30 

99.60 

99.50 

99.40 

99.50 

Vision  Gauge 

11 

. 014 

. -026 

. 0.04 

. -026 

. -026 

. 024 

. -016 

044 

. 024 

. -026 

. 024 

Angle 

96.62 

96.60 

96.40 

96.80 

96.80 

96.40 

96.30 

96.60 

96.40 

96.80 

96.30 

96.20 

Vision  Gauge 

12 

. -0.02 

. -022 

. 018 

. 018 

. -022 

. -022 

-6.62 

. J 022 

. 018 

. -022 

. -042 

833.00 

832.40 

832.60 

832.10 

832.15 

832.30 

832.55 

832.10 

832.35 

832.90 

832.45 

832.25 

Microscope 

-1 .00 

. -060 

. -040 

. -0.90 

. -025 

. -070 

. -045 

. -020 

. -025 

. -010 

. -055 

. -075 

102.00 

102.00 

101 .85 

102.05 

101 .90 

101 .75 

102.30 

102.20 

102.30 

102.15 

102.45 

102.10 

Microscope 

+/-2.00 

0.00 

-0.15 

0.05 

-0.10 

-0.25 

0.30 

0.20 

0.30 

0.15 

0.45 

0.10 

144.00 

143.10 

143.05 

143.00 

143.20 

143.15 

143.05 

143.20 

143.25 

143.05 

143.50 

143.20 

Microscope 

+/-1 .00 

-0.90 

-0.95 

-1 .00 

-0.80 

-0.85 

-0.95 

-0.80 

-0.75 

-0.95 

-0.50 

-0.80 

Mat.  thickness 

29.53 

29.53 

29.61 

29.53 

29.53 

29.53 

29.53 

29.92 

29.53 

29.53 

30.51 

Maq 

x32 

x32 

x32 

x32 

x32 

x32 

x32 

x32 

x32 

x32 

x32 
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—♦—#4  (+0.40:  0.05  -  0.40  mil) 
—•—#5  (+7-0.5:  -0.03  -  0.40  mil) 
— #6  (+-/-0.5:  -0.05  -  0.24  mil) 
-■-#13  (-1.0:  -0.9- -0.1  mil) 


Measured  Dimensions 
Holes 


Edge 

10.00 

9.84 

9.88 

10.00 

10.75 

10.12 

9.50 

10.00 

10.12 

10.31 

10.35 

9.96 

Top 

6.00 

6.33 

6.20 

6.20 

6.13 

6.16 

6.36 

6.47 

6.49 

6.46 

6.12 

6.08 

Dia 

+/-0.50 

0.33 

0.20 

0.20 

0.13 

0.16 

0.36 

0.47 

0.49 

0.46 

0.12 

0.08 

14.80 

15.08 

14.96 

14.65 

14.53 

14.65 

14.84 

14.96 

14.88 

14.80 

15.08 

ym 

14.76 

14.80 

14.69 

14.80 

14.53 

14.92 

14.80 

15.00 

14.45 

14.92 

15.35 

Bott 

6.00 

6.26 

6.36 

6.29 

6.27 

6.13 

6.44 

6.34 

6.45 

6.49 

6.08 

6.11 

Dia 

+/-0.50 

0.26 

0.36 

0.29 

0.27 

0.13 

0.44 

0.34 

0.45 

0.49 

0.08 

0.11 

14.96 

15.24 

14.76 

15.24 

14.57 

14.76 

14.96 

15.04 

15.04 

14.65 

14.84 

ym 

14.57 

14.57 

14.88 

14.65 

15.00 

14.96 

14.65 

14.96 

14.21 

14.88 

15.51 

B.C. 

135.00 

134.84 

135.15 

135.00 

135.20 

135.05 

135.10 

135.05 

135.15 

135.10 

135.05 

135.15 

Microscope 
V.  Gauge  64x 

Microscope 
Microscope 
V.  Gauge  64x 

Microscope 

Microscope 

Microscope 


top  hole  diam.  (+7-0.5:  0.08  -  0.49  mil) 
bottom  hole  diam.  (+7-0.5:  0.08  -  0.49  mil) 
from  top  (+7-1 :  -0.5  -  0.75  mil) 
between  holes  (+7-0.5:  -0.16  -  0.20  mil) 
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Measured  Dimensions 
Test  Bench 


I  Part# 
Dim  j  Req 

. 18 . 

Mill  3 . 

. 19 . 

. Miii'3 . 

155.00 

154.90 

154.70 

+/-1  .00 

-0.10 

-0.30 

„  30.00 

30.30 

30.45 

r  +20Q 

0.30 

0.45 

3  „ . 1000 

10.35 

10.00 

3  r  +0.40 

0.35 

0.00 

4***  „ . 11,00 

11.20 

11.00 

4  r  +0.40 

0.20 

0.00 

5  L  5.00 

5.35 

5.15 

Rad  +/-0.5 

0.35 

0.15 

6  R  5.00 

5.08 

5.00 

Rad  +/-0.5 

0.08 

0.00 

,  30.00 

29.45 

29.15 

+/-5.00 

-0.55 

-0.85 

65.00 

64.35 

64.45 

8  +/-5.00 

-0.65 

-0.55 

105.00 

105.80 

106.35 

+/-5.00 

0.80 

1.35 

90  00 

89.15 

89.45 

10  . <7.5.00 

-0.85 

-0.55 

Angle  99.46 

99.60 

99.50 

11 

0.14 

0.04 

Angle  96.62 

96.30 

96.40 

12 

-0.32 

-0.22 

833.00 

832.50 

833.00 

13  :t:oo 

-0.50 

0.00 

16  . 9.4,00. 

93.35 

93.50 

-2.00 

-0.65 

-0.50 

h  (FW050901 01  _test_bench) 
Method 


Microscope 

Microscope 

Microscope 

Vision  Gauge 

Vision  Gauge 

Microscope 

Microscope 

Microscope 

Microscope 

Vision  Gauge 

Vision  Gauge 

Microscope 

Microscope 


Part# 

. 18 . 

. 19 . 

Dim  j 

Req 

. Mill  3 . 

. iviiii  3 . 

16 

94.00 

93.35 

93.50 

-2.00 

-0.65 

-0.50 

255.00 

255.25 

255.20 

+/-5"00 

. 6.25 

. CL20 

101 .00 

100.05 

100.40 

+/-1 .00 

-0.95 

-0.60 

19  • 

59.00 

58.95 

59.00 

+/-5.00 

-0.05 

0.00 

20  • 

32.00 

31 .35 

31 .50 

+/-5.00 

-0.65 

-0.50 

21  • 

19.00 

19.70 

19.70 

+/-5.00 

0.70 

0.70 

Mat.  thickness 

30.51 

30.51 

Mag 

x32 

x32 

Method 

Microscope 

Microscope 

Microscope 

Microscope 

Microscope 

Microscope 


Dim 

Part# 

Req 

"ZKZ 

ZZZZ 

Edge 

10.00 

10.24 

10.24 

6.00 

6.13 

6.15 

Top  DI9  • 

™;ao5o 

. 6.13 

0.15 

15.35 

15.24 

ym 

. 1 5.28 

. 1524 

6.00 

6.20 

6.10 

□Ott  UI8  • 

+/-0'50 

. 616 

0.10 

Sym 

15.08 

15.24 

15.16 

15.24 

B.C. 

135.00 

135.25 

135.25 

Microscope 
Vision  Gauge 

Microscope 
Microscope 
Vision  Gauge  <E 

Microscope 

Microscope 

Microscope 
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part  #1  part  #10  part  #15 


xlOOO 


6.1-13 


part  #1 


part  #10  part  #15 


4-5  parts  per  tool 


6.1-14 


6.1-15 
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7.  Control  Synthesis 


7.1  Position  Control  of  SMA  Actuation 

presented  by  N.  Lechevin,  DRDC-Valcartier 

A  practical  digital  controller  was  designed  with  the  objective  of  controlling  the  effector 
position  of  an  SMA  actuation  mechanism  over  a  range  of  1  mm  and  within  a  bandwidth  of  [0  Hz, 
1  Hz].  The  plant  under  digital  control  comprised  the  compliant  flow  effector,  the  SMA  wires,  the 
power  drive  and  the  electronics.  To  meet  the  stringent  positioning  and  rate  requirements,  a  two- 
step  variable  structure  set-point  regulator  was  proposed.  The  controller  is  comprised  of  two  parts. 
First,  to  guarantee  fast  transients,  a  bang-bang  controller  is  activated  whenever  the  regulation 
error  is  above  a  fixed  threshold.  Second,  a  linear  digital  controller  is  switched  on  when  the  error 
reaches  the  threshold.  The  linear  digital  controller  provides  smooth  convergence  of  the  closed- 
loop  system  to  the  steady  state  and  prevents  high  frequency  chattering,  which  is  typical  of  bang- 
bang  control. 

Since  the  available  micromechanical  model  of  the  SMA  actuator  was  not  intended  for 
control  studies,  a  discrete-time  parameter  identification  of  the  actuated  flow  effector  was  carried 
out  prior  to  control  synthesis.  With  the  discrete-time  linear  plant  model  obtained,  the  design  of  a 
digital  controller  followed.  The  closed-loop  discrete-time  control  system  had  to  satisfy  1)  closed- 
loop  pole-placement  and  2)  steady-state  error  requirements.  Experimental  results  showed  that  the 
two-step  variable  structure  controller  satisfies  the  required  specifications,  providing  response 
times  below  0.3  sec  with  overshoot  of  about  5%  of  the  steady  state  value. 
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1-  Control  System  for  SMA-based  Effector  Positioning 


1-  Control  System  for  Effector  Positioning 
(Cont’d) 


LVDT:  linear  variable  displacement  transducer; 
SMA  :  shape  memory  alloy. 
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2.  Objectives 

>  Robustness  with  respect  to  parametric  uncertainties 
such  as  friction  and  the  aerodynamic  force. 

>  Compensate  for  the  hysteretic  behavior  of  the  actuator; 

>  Good  tracking  performance  within  [0  Hz,  1  Hz]; 

>  Zero  steady-state  error  =>  avoid  chattering; 

>  Displacement  Range  :  1  mm 


3.  Previous  Studies  in  Open  Literature 


Hysteresis  model-based  Control:  Preisach  model  and  the  like 
to  compensate  for  hysteresic  nonlinearity  -  [0,0.1  Hz] 


■  Adaptive  control  :  [Webb  et  al.,  98],  [Kokotovic  et  al.,  96] 

■  Stochastic  approach  +  PID  controller  +  PWM:  [Majima  etal,  01]; 


Non-model-based  Control 


■  Dissipativity-based  PD  controller:  [Gorbet  etal.,  97-01] 

■  Time-delay  control:  [Lee  et  al.,  04  ]; 

■  Quasi-linearization  by  open-loop  compensation  with  phaser: 
[Cruz-Hernandez  and  Hayward.,  01,  05]; 


■Variable  structure  control:  [Grant  and  Hayward.,  97]; 

o  Sensitivity  to  model  uncertainties; 
o  Mechanical  load  is  not  considered,  [Gorbet,  01], 


Drawbacks: 
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4-  Variable  Structure  control 


4-  Variable  Structure  Control 

The  design  of  the  linear  DT  control  law  is  based  on  a  linear 
approximation  with  a  pure  delay  of  the  SMA  actuator. 

■  Quasilinearization  by  open-loop  compensation  with 
phaser:  [Cruz-Hernandez  and  Hayward.,  01,  05] 

=>  nonparametric  identification:  estimation  of  the 
phase  lag  in  the  Bode  diagram; 

=>  design  of  a  phaser  to  compensate  for  the  delay; 

■  1.  Parameter  identification  (ARX  model) 

2.  A  LTI  DT  controller  is  designed  from  1. 
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5-  Parameter  Identification 


5^_  Parameter  Identification  (Cont’d) 

a)  Experiment 

INPUT  ;  Band-limited  Chirp  Signal  v  (volt)  0.3  Hz  <  /  <1  Hz 
>  Each  test  is  performed  with  a  constant  amplitude:  Va 
OUTPUT:  Measured  displacement  (dj  stored  in  the  board 
INPUT  and  OUTPUT  are  stored  for  each  Va 
Va  is  swept  from  1  V  to  4.5  V 


b)  Off-line  parameter 
identification 


FOR  each  Va  GET  (v,  dj 
FOR  na=  1:5 

FOR  na=  1:5 

FOR  d=l:15 
(ap...,ana,b 
END  END  END  END 


ARX  model:  y{t)  +  aly{t-\)  +  ...  +  amy(t-na)  = 


y-dm, 

u:  v  , 


hu(t  -  nk )  + ...  +  b bu(t  -d-nb  + 1) 

d:  pure  delay 

na:  number  of poles  ,  nb-1  :  number  of  zeros 


p..., bnb )  =ARX(v,  dm, na,nb,d) 


KEEP  coefficients 
which  provide 
a  fitting  ratio  >  80% 
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6-  DT  Control 


Output  of  the 

outer-loop  Aero- 

Roll  control  Antagonistic  actuation  dynamics 


Closed-loop  system 


DT  Control  law  objectives: 

□  Fast  and  Accurate  set-point  regulation 

□  Achieve  smooth  transition 

►  Pole  placement  of  the  closed  loop 

►  Low-pass  filter 

►  Update  initial  condition  of 
controller  integrators 


6-  DT  Control  (Cont’d) 

Output  of  the 

outer-loop  Control  law  with  SMA  model 

77  7  an  integral  action 

roll  control 


COMPUTATION  of  the  Control  Law  matrices:  direct  digital  method 

[Astrom  and  Wittenmark,  90] 

□  Include  an  integrator  in  front  of  the  plant; 

□  Closed-loop  pole  placement ; 

□  Solve  a  Diophantine  equation  system  to  obtain 
coefficients  of  polynomials  P(q~r)  and  L(q-])  . 
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7-  Experimental  Results 

*♦♦  Sample  period  :  5  ms  (time  constants  of  the  heat  transfer  dynamics  and 

of  the  compliant  mechanism  dynamics  are  larger 
than  100  ms); 

❖  deg  B  =  1  ,  deg  A  =  2; 

❖  Pole  placement:  Sj  =  s2  =  -3  ,  s3  =  s4  =  s5  =  -30; 

v(k)  _  l-1.5676q  1  +  0.5918q 2 

~e(k)~  3.4062 -6.5593q  1  +3.1582q'2 
s  -  0.15  mm 
V  -  3  volts 

♦>  PWM  voltage  source  control:  fsw  =  400  Hz;  0V  <  v  <  3; 

❖  SMA  Wires:  length  =  5  in,  diameter  =  0.004  in; 

♦♦♦Set-point  command:  0  mm  <d*<  1  mm 

f  =  0.25  Hz  and  0.5  Hz. 
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8-  Conclusion 

>  Variable  Structure  Controller  is  devised  to  achieve 
robust  control  and  fast  positioning  of  the  effector; 

>  A  linear  Discrete-time  law  based  on  model  identification 

is  designed  to  achieve  zero  steady-state  error  without  chattering. 

Future  Work:  increase  actuator  response  times 

□  Investigate  cooling  phase:  smaller  diameter  wires 
=>  tradeoff  between  diameter  and  mechanical  load; 

□  Increase  the  moment  arm  of  the  bottom  wire. 
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7.2  Nonlinear  Control  Methods  for  Delta  Wing  Systems  with  SMA 
Actuator  Dynamics 

presented  by  B.  Gordon,  Concordia  U. 

This  presentation  investigates  the  application  of  various  nonlinear  control  methods  for 
delta  wing  systems  using  shape  memory  alloy  (SMA)  actuator  models  recently  developed  at 
DRDC.  Radial  basis  function  (RBF)  neural  network  control,  hybrid  RBF  control,  receding 
horizon  control  (RHC),  and  linear  parameter- varying  (LPV)  control  approaches  are  applied  to  the 
system.  The  performance  and  robustness  of  the  proposed  control  methods  in  the  presence  of 
uncertainty  and  SMA  actuator  time  delays  are  examined.  Furthermore,  the  applicability  of  the 
techniques  to  future  experimental  work  and  possible  implementation  issues  are  discussed. 
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Nonlinear  Control  Methods 
for  Delta  Wing  Systems  with 
SMA  Actuator  Dynamics 


Mehrdad  Pakmehr,  Ming  Yang,  Hojjat  Izadi,  and  B.  W.  Gordon 
CIS  Laboratory,  Concordia  University 
April  2006 
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Vortex  Coupled  Delta  Wing  System 


■  Aerodynamic  loads  are  dependent  on  the  vortex 
breakdown  location  over  the  delta  wing 

■  Left  and  right  primary  vortices  need  to  be  considered 

Xvbl  =  Xsl(0(t))  +  Xq(<j>(t))  ■  Xsl(<f)(t))  +  f  Xu(t  -  T)<i>{T)dT 

Xvbr  =  Xsr(</>(t))  +  X^(t))  ■  Xr{<t>{t))  -  j'  Xu(t  ~  T^>(T)dT  3 


Delta  Wing  State  Space  Model 

■  State  space  model  (Liu  and  Gordon,  2004) 

x\(t)  -c-xl(t) 

xi  (t)  =  -c  ■  xi  (t) + x4  (t)  +(x4 

X3(t)  =  X4  (() 

X4(t)  =  -Cl(Xvbl,  Xvbr)  ■  Q  -  fc  ■  X4 (t)  +  1 II-  Ui(t) 
y(t)  =  X3(0 

Xvbl  =  Xsl(x3)  + Xsl(x3)-kq(t)- X4(t)  + Cl(t)- Xl(t)  +  Ul(t) 
Xvbr  =  Xsr(x 3)  +  Xsr(x 3)  •  kq(t)  •  X4 (t)  ~  a{t )  •  Xl(^)  +  Ul{t) 

■  4th  order  nonlinear  model  with  time  delay  T 
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Objectives 


■  To  develop  nonlinear  tracking  controllers  for  vortex  coupled 
delta  wing  dynamics 

■  Controller  must  compensate  for  unknown  parameters, 
uncertainty,  and  time  delays 

■  Combine  the  nonlinear  control  methods  with  shape  memory 
alloy  (SMA)  actuator  dynamics 


Multiple  Input  Controller  Structure 

u3(t )  =  u(t ) 

«i  (0  =  P.bMt) 

“2(0  =  -PvbAt) 

■  u(t)  is  the  main  control  input  which  is  the  same  as 
what  has  been  presented  in  the  previous  section 

■  The  other  two  control  inputs  have  been  defined  as 
perturbations  in  the  left  and  right  vortex  breakdown 
positions 

■  The  relative  weighting  of  the  vortex  inputs  is  given 
by  Pvbi 
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Stabilizing  Adaptive  Control  Approach 

A  stabilizing  adaptive  controller  was  developed  based  on 
existing  works  (Foda,  1998) 

The  controller  is  based  on  output  measurements 
uncertainty  bounds  are  assumed 


i(r)  =  A0x(t)  +  D0x(t  - r/)  +  B0[udp(xt ,t)  +  h(xt , Xvb,t)] 

Up(xt’t)  =  -J“dHoy(t) 


Bd  =  SBd  +FoJKO| 
x.  (6)  =  v{0)  9  e  [-77 ,0] 


9d(t0)  =  9d 


Output  feedback  Matrix  H0 

1-  In  our  case  as  the  rows  of  c„  are  linearly  independent  so  (c0c0r) 1 
does  exist,  and  hence  (c0c0+)=i,  so  we  can  obtain  h0  as  follows: 

h0=b0tp0c ; 

We  can  obtain  H0  only  if  the  matching  condition  has  a  feasible 
solution. 

2-  Defining  a  constrained  linear  matrix  inequality  (LMI)  as  follows: 


A>  0  <  0 


subject  to  constraints 


H0C0=B0‘P0 


Hn  <Ha<Hn 

U  min  U  U  max 


this  can  be  solved  using  convex  optimization  tools,  like  YALMIP 
Toolbox,  Matlab  LMI  Toolbox  or  LMITools,  operating  in  Matlab 
environment. 


Flowchart  of  practical 
application  of  delta 
wing  stabilization 
control  process 


Schematic  of  the  adaptive  stabilizing  controller 


To  improve  the  control  task  we  define  the  following 
mechanisms  for  adaptation  gain: 

j'M/  <  t  &  fa  <  0  =>//„  = 

[if  >  Mmax  &Aa>0=>Md=  /'max 
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Modelling  Uncertainty  Testing 

To  test  the  robustness  of  the  method  we  defined  the 
following  disturbance  in  c,: 

SC,  =  Kean  sin(®0 

where  the  Amec,n  is  a  mean  value  of  a  is 

A<a  +  Z>||x*||,  A  =  \ct  -c,\ ,  x]  =x/xsi 

a  =  0.00047,  b  =  0.00258 
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Simulation  Results  for  Stabilizing  Controller 


output  time  history 
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control  input  time  history 


u  (control  input)  time  history 


variation  of  the  controller  gain  and  adaptation  law 
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Left  and  Right  vortex  breakdown  time  history 


Variation  of  the  controller  gain 
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Hybrid  Adaptive  Neural  Network  Control 


This  approach  is  proposed  to  address  the  multiple  equilibrium  problem 
and  estimation  of  Cl  from  multiple  experimental  data  sets 
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Hysteresis  switching  logic  for  hybrid  control 


021  —  01  K  012  —  0]  +  hs  032  02  As  023  02  As 

— ■ — ^ — ■ — — Ml — ■ — 

01  02 


■  Hysteresis  switching  logic  is  proposed  to  reduce  chattering 
and  limit  cycling  between  the  finite  states  of  the  hybrid  system 


Simulation  Results  for  Hybrid  Controller 

A  =  -20  [deg], (4  =  20  [deg],  hs  =  4[deg] 


Sliding  time  history  Error  time  history 
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Neural  Network  Adaptive  Control  Combined  with 
SMA  Model  Developed  at  DRDC 

Vortex  coupled  delta  wing  model: 

Xx{t)=CX2{t) 

x2  ( t )  =-cxl  ( t )  -  sd  x2  (t) + x4  (t) +x4(t- T) 

x3  (t)  =  -sd  x3  (t)  +  x4  (t) 

x4(t)  =  -C,(K)Q-fc  x4(t) +hlul(t)/Iw 

K  =  [Xvbl(t),Xvbr(t)] 

XM  ( t )  =  Xsl  (x(tj) + xsl  (40)  kq  (0  4  (0 + 40  4  (0 + ysma,  (0 
XvJt)  =  XJx(t)) + XJx(t))kq  (?)  x4(t)  -  a(0  X,  (?) + ysmar(t) 

y(t)  =  xs(0 

19 


Simulink  block  diagram  of  continuous  time 
approximation  of  SMA  control  loop 


Output  of  the 
outer-loop 
Roll  control 


CT  Control 

Antagonistic  actuation 


Aero¬ 

dynamics 


t 


j 


*Lechevin,  N.,  Boissoneault,  O.,  and  Rabbath,  C.  A.,  “Control  of  Antagonistic  Actuator  by  Means  of 
Shape  Memory  Alloy  -  Identification-Based  Control”,  Internal  Report,  Defence  Research  and 
Development  Canada  (DRDC),  Valcartier,  Quebec,  Canada,  January  5,  2006. 
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Continuous  SMA  Control  System  Model 

Low  pass  filter 

V  0)  =  AwxiPf(  0  +  ^/«,(0 
j/P/(0  =  Cip/Xip/iO + DiPfui(t ) 

II 

Control  law  with  an  integrator 

(0  =  +  BcontUconM ) 

+  Dcon,ucont{t ) 

Ucont(t)  =  (ylpf(t)-ysma(t )) 

Identified  SMA  model 

i:  (t)  =  A  x  (t)  +  B  u  (t) 

V  (t)  =  C  X  (t)  +  D  U  ( t ) 

sma  V  /  sma  sma  V  /  sma  sma  V  / 

“„™(0  =  ^OwCO) 
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Simulation  of  the  combined  control  system 


Roll  angle  time  history 
(plant  output  and  the  command  signal) 


Error  time  history 


Left  and  Right  SMA  control 
loop  output  time  history 


Overall  control  input  time  history 


23 


Linear  Parameter  Varying  Sliding  Mode 
Control  (LPVSMC)  Approach 

LPVSMC  drives  the  state  trajectory  of  a  state  delayed  system  onto 
a  parameter-varying  sliding  surface: 

s(o(t))  =  H  (a(t  ))  z  w  +  L  (o(t))  r  -  L  (o(t))  C  (a(t))  z,+z2 


LPVSMC  input:  u(o(t))  =  |u-(a(t))_u"(a(t))’  s(a(t,)  '  6(a(t»>0 

[ueq(o(t))+u„(o(t)),  s(a(t))TB(a(t))<0 

u  n(o(t»  =  7 - 11  \  | - 

s(o(t))TB(o(t))  +  £  24 
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LPVSMC  Approach 


After  the  sliding  surface  is  reached  the  state  asymptotically 
approaches  the  desired  trajectory 


The  stability  of  the  state  delayed  system  on  the  sliding  manifold  can  be 
proven  using  a  parameter  dependant  Lyapunov-KrasovskiT  functional 


V  =  V1(o(t))  +  V2(o(t))  +  V,(o(t)) 
V1(o(t),t)  =  zT(t)P(o(t))z(t) 

v2  t)  =  Jt'_h  z T  (e)R(o(e))z(e)de 

V3  (o(t),  t)  =  jt'  h  j 'z  T  (e)Q(c(e))z(0)d0ds 
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Simulation  of  LPVSMC  with  1st  order  +  delay 
approximation  of  SMA  actuator  dynamics 


The  LPVSMC  approach  is  robust  to  model  uncertainty  and  time  delays 

The  performance  is  much  better  than  regular  LPV  control  for  the 
vortex  control  problem 
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Receding  Horizon  Control  Approach 

nominal  model:  x(t)  =  f(x(t),u(t))  x(0)  =  x0 

t+T 

RHC  optimization  problem:  min  jq(x(T),u(r))dr  +  V(x(t  +  T)) 

t 

Subject  to  input  or  state  constraints:  lbt  <  <  ubt 

Closed  loop  RHC  solution:  x(r)  =  f{x(z),u\z)) 

t  <=[t,t  +  S),  0  <5  <T 

RHC  control  is  well  suited  to  achieve  high  performance  with  systems 
subject  to  input  saturation  and  other  operating  constraints 
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Simulation  of  RHC  with  1st  order  +  delay 
approximation  of  SMA  actuator  dynamics 


xl  ■  xl 


■r - ‘ ‘ ^ ^ - 1 

0  1  2  3  4  5  6  J _ ,  ,  ,  , _ I 

Time,  sec  0  1  2  3  4  5  6 

Time,  sec 


Simulation  studies  indicate  RHC  control  is  robust  to  model  uncertainty 
and  time  delays  for  the  vortex  control  problem 

However,  no  theoretical  guarantees  are  currently  available  for  this  class 
of  problems  28 
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Comparison  of  Nonlinear  Control  Approaches 


Tracking 

performance 

Saturation 

effects 

Chattering 

effects 

good 

high 

none 

excellent 

medium 

high 

excellent 

medium 

high 

good 

low 

medium 

excellent 

none 

none 

LMI 

Adaptive 

Stabilization 


Neural  network 
adaptive 
control 


Hybrid  neural 
network 
adaptive 
control 


Linear 
parameter 
varying  sliding 
mode  control 


Receding  excellent 
horizon  control 


*  Observed  from  simulations,  but  no  theoretical  guarantees  currently  exist  29 


Conclusions  and  Future  Work 

■  Nonlinear  tracking  controllers  for  vortex  coupled  delta  wing 
dynamics  were  developed  and  tested  on  simulations 

■  The  nonlinear  control  methods  were  combined  with  shape 
memory  alloy  (SMA)  actuator  dynamics 

Possible  Future  Work: 

■  Full  integration  of  the  proposed  vortex  control  methods  with  the 
SMA  actuator  control  approach 

■  Experimental  implementation  and  investigation 

■  CFD  investigation 
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